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INTRODUCTION 

The  purpose  of  this  research  project  was  to  investigate  the  effect  of  small  scale 
fluid  mechanics  on  microbial  uptake  of  low  and  high  molecular  weight  dissolved  organic 
compounds.  Our  hypothesis  was  that  microbial  :>ptake  of  large  molecular  weight 
compounds,  under  specific  conditions,  could  be  increased  by  fluid  shear  and  advective 
flow  past  the  microorganisms.  This  hypothesis  was  based  on  a  biohydrodynamical 
uptake  model  derived  from  mass  transfer  and  coagulation  theories  (Logan  and  Hunt 
1987,1988).  From  equations  used  to  describe  mass  transfer  to  microbes,  and  coagulation 
of  microbes  with  substrate,  it  was  predicted  that  substrate  uptake  could  be  increased  by 
microorganisms  in  high  shear  environments  versus  microorganisms  in  nearly  stagnant 
flow.  However,  the  relative  advantages  of  growth  within  high  shear  environments  was 
calculated  to  increase  with  increasing  molecular  weights  of  organic  matter.  Thus,  for 
low  molecular  weight  compounds,  increased  uptake  with  fluid  shear  should  be 
undetectable.  The  majority  of  uptake  studies  reported  in  the  marine  literature  have 
focused  on  easily  degraded  compounds  with  small  molecular  weights,  such  as  glucose  and 
amino  acids,  despite  research  which  indicates  large  molecular  weight  organics  can 
account  for  a  substantial  fraction  of  dissolved  organics.  Therefore,  researchers  would 
not  have  observed  an  effect  of  fluid  shear  on  uptake. 

Experiments  were  proposed  to  study  microbial  uptake  in  two  different  fluid 
environments:  suspended  in  fluid  shear  and  fixed  in  a  uniform  flow  field.  The  plan  of 
study  included:  (a)  defining  pure  culture  microbial  uptake  as  a  function  of  substrate  size 
and  fluid  environment  using  radiolabeled  substrates  to  validate  the  biohydrodynamic 


model;  (b)  continued  development  of  the  biohydrodynamic  model  and  mass  transfer 
correlations;  (c)  development  and  application  of  techniques  to  measure  uptake  of  large 
molecular  weight  compounds,  such  as  protein.  The  results  of  this  research  is 
summarized  below.  Details  of  our  results  are  contained  in  the  attached  research 
publications. 
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SUMMARIZED  RESEARCH  RESULTS 


INCREASED  MASS  TRANSFER  TO  MICROORGANISMS  WITH  FLUID  MOTION 

The  effect  of  fluid  flow  and  laminar  shear  on  bacterial  uptake  was  examined  in 
the  laboratory  under  a  wider  range  of  conditions  than  could  be  experienced  by 
unattached  and  attached  marine  bacteria.  During  the  first  phase  of  research,  we 
examined  leucine  uptake  by  pure  cultures  of  suspended  (unattached)  bacteria  (Z. 
ramigera).  Since  our  mathematical  analysis  assumed  bacterial  transport  systems  were  not 
saturated,  we  first  examined  the  kinetics  of  leucine  uptake  by  this  bacterium.  We 
measured  leucine  uptake  at  three  difference  incubation  times,  at  concentrations  between 
0.5  /ig  l"1  to  100  mg  I"1.  Leucine  uptake  was  linear  with  concentration  over  a  range 
spanning  6  orders-of-magnitude.  As  incubation  times  increased,  the  rate  of  uptake 
decreased  suggesting  feedback  from  an  intracellular  leucine  pool.  For  incubation  times 
of  0.5,  1  and  7  minutes,  the  slopes  (±SD)  of  uptake  versus  concentration  (on  log-log 
axes)  were  7.0+0.02,  6.4±0.02,  and  5.2+0.03  x  lO^min-1,  respectively.  The  r2  values  for 
all  lines  were  greater  than  0.98.  These  results  indicate  that  saturation  kinetics  were  not 
exhibited  by  2L  ramigera  over  a  wide  range  of  leucine  concentrations. 

Several  flow  experiments  conducted  with  Z.  ramigera  held  on  0.2  fi m 
polycarbonate  filters  indicated  leucine  uptake  by  cells  fixed  in  a  flow  field  was  55-65% 
greater  than  leucine  uptake  by  suspended  cells.  In  these  experiments,  leucine  uptake 
increased  with  fluid  velocity  at  velocities  less  than  ~1  mm  s“*.  Little  or  no  increase  in 
leucine  uptake  was  observed  at  higher  velocities  between  1-2.5  mm  s  .  On  occasion, 
some  cells  were  completely  inactivated  at  higher  flow  velocities  and  did  not  take  up  any 
leucine,  probably  due  to  cell  disruption  at  vacuums  of  600  mm  Hg  needed  to  obtain 
higher  velocities.  Fluid  velocities  below  1  mm  s~^  were  obtained  at  vacuums  below  250 
mm  Hg,  a  range  indicated  by  Fuhrman  and  Bell  (1985)  to  prevent  significant  release  of 
amino  acids  from  marine  microorganisms. 
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Fluid  flow  experiments  using  0.45  cellulose  acetate  and  0.2  /tm  aluminum  oxide 
filters  to  hold  bacteria  were  characterized  by  large  variations  in  bacterial  uptake  about 
mean  values.  Cellulose  acetate  filters  retained  large  concentrations  of  [^HJleucine  label 
during  flow  experiments  which  could  not  be  removed  by  filter  rinsing.  This  substantially 
increased  the  amount  of  abiotic  adsorption  of  leucine  to  filters.  These  increased 
background  leucine  concentrations  produced  large  variations  in  uptake  for  a  given  fluid 
velocity.  Based  on  average  uptake,  the  effect  of  fluid  motion  on  uptake  was  similar  to 
that  observed  with  the  polycarbonate  (surface)  filters.  A  maximum  average  increase  of 
60%  was  consistent  with  previous  results  using  the  polycarbonate  filters.  However,  the 
larger  errors  associated  each  measurement  did  not  result  in  a  significant  correlation  of 
fluid  velocity  with  uptake. 

The  effect  of  fluid  motion,  predicted  from  a  mass  transfer  analysis,  was 
compared  with  experimental  results.  Mass  transfer  calculations  were  based  on  the 
assumptions  that  leucine  concentrations  at  the  cell  surface  were  much  less  than  the  bulk 
concentration,  that  uptake  was  not  zero-order  (saturated),  and  that  flow  past  cells  was 
completely  uniform.  At  a  flow  velocity  of  1  mm  s  ,  three  different  mass  transfer 
correlations  predicted  increases  22,  65  and  110%,  versus  the  observed  increase  of 
55-65%.  The  correlation  used  by  Munk  and  Riley  (1952)  consistently  overestimated  the 
effect  of  fluid  motion.  The  mass  transfer  correlation  developed  by  Logan  and  Alldredge 
(1989)  gave  the  best-fit  to  the  experimental  results  of  this  study. 

The  effect  of  fluid  shear  on  leucine  uptake  by  Zu  ramieera  was  examined  by 
suspending  cells  in  a  laminar  shear  field  generated  using  a  rotating  cylinder  device 
constructed  for  this  research.  Within  experimental  error,  fluid  shear  did  not  increase 
leucine  uptake  by  bacteria.  Moreover,  there  was  a  consistent  trend  of  decreasing  uptake 
with  shear  that  suggested  that  higher  laminar  shear  rates  impaired  leucine  uptake  by 
suspended  bacteria.  This  result  confirmed  the  prediction  of  the  mass  transfer  analysis 
indicates  that  fluid  shear  should  not  increase  leucine  uptake  by  more  than  5%  at  shear 
rates  below  50  s" 1 . 


-4- 


Conclusions:  From  both  experimental  results  and  theoretical  considerations,  fluid 
shear  is  not  expected  to  substantially  increase  microbial  uptake  by  unattached  bacteria. 
Uptake  experiments  with  2^  ramieera  showed  that  fluid  shear,  and  therefore,  fluid 
mixing  intensity,  did  not  increase  leucine  uptake  by  bacteria  suspended  in  a  laminar 
shear  field.  Fluid  flow  past  bacteria,  however,  can  substantially  increase  bacterial  uptake. 
Leucine  uptake  by  Z L  ramieera  increased  by  an  average  of  55-65%  at  flow  rates  up  to  1 
mm  s'1.  In  order  for  cells  to  be  present  in  this  type  of  flow  field  under  natural 
conditions,  they  must  either  be  attached  to  a  surface  or  attached  to  other  cells  in 
permeable  aggregates.  We  believe  the  effect  of  fluid  motion  on  bacterial  kinetics  is 
attributable  to  compression  of  the  concentration  boundary  layer  surrounding  surface 
transport  proteins.  The  rate  of  leucine  uptake  suggests  that  binding  proteins  used  to 
transport  leucine  into  the  cell  can  occupy  approximately  1%  of  the  cell  surface  area.  A 
detailed  description  of  this  research  is  contained  in  Logan  and  Dettmer  (1990). 

INCREASED  UPTAKE  OF  DISSOLVED  ORGANICS  BY  MARINE  BACTERIA  AS  A 
FUNCTION  OF  FLUID  MOTION 

The  purpose  of  this  investigation  was  to  examine  if  fluid  motion  affected  uptake 
by  natural  assemblages  of  marine  bacteria  in  the  same  manner  as  a  pure  culture.  Fluid 
motion  should  only  increase  uptake  when  uptake  is  not  already  saturated  and  at  the 
maximum  rate.  We  measured  the  uptake  of  two  model  compounds  (leucine  and  glucose) 
by  marine  bacteria  fixed  in  a  uniform  flow  field  and  suspended  in  a  laminar  shear  field. 

All  samples  were  collected  from  surface  waters  in  the  Roosevelt  Inlet  on  the 
Delaware  Bay  (Lewes,  DE)  during  August,  1989  through  the  assistance  of  David  L. 
Kirchman  (University  of  Delaware,  Lewes).  Surface  samples  were  gravity-filtered 
through  0.8  nm  Nucleopore  filters  to  remove  bacterivores  and  autotrophs.  Samples  were 
either  immediately  used  in  experiments  or  kept  overnight  (aged)  in  the  dark  at  room 
temperature.  Leucine  and  glucose  uptake  by  bacteria  was  determined  using 
4,5-[-*H]leucine  and  [^HJglucose. 
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Fluid  motion  was  predicted  to  increase  microbial  uptake  only  if  transport  were 
not  saturated.  During  our  initial  experiments,  leucine  transport  was  probably  saturated 
in  surface  waters  of  Roosevelt  Inlet  since  the  leucine  concentration  was  45  nM  (August 
7)  and  Ks  for  leucine  uptake  is  probably  <1  nM.  To  reduce  leucine  concentrations  and 
obtain  non-saturated  uptake  kinetics,  we  aged  the  bacterial  size  fraction  (<0.8  ^m 
filtrate)  for  24  hr  in  the  dark;  in  th.e  absence  of  leucine  input  from  other 
microorganisms,  bacterial  uptake  reduced  leucine  concentrations  to  1.7  nM.  Leucine 
concentrations  in  particle-free  seawater  were  near  the  detection  limit  (<1  nM),  indicating 
that  the  filtered  water  used  for  the  flow  experiments  was  not  contaminated. 

At  low  leucine  concentrations  (aged  samples),  we  found  that  fluid  flow  past  cells 
increased  [^H]leucine  uptake  by  a  factor  of  8  compared  to  uptake  at  low  fluid  velocities. 
[^H]Leucine  uptake  increased  from  0.65  pg  ml-*  min-*  at  3  m  d'*  to  a  maximum  of  5.4 
pg  ml-*  min-*  at  -20  m  d-*.  At  higher  fluid  velocities,  [^H]leucine  uptake  decreased 
to  an  average  of  4.5  pg  ml-*  min-*  for  fluid  velocities  between  20-70  m  d-*.  The 
maximum  rate  of  leucine  uptake  in  flow  experiments  was  63%  less  than  the  average 
uptake  of  suspended  cells  (14.5±0.7  pg  ml-*  min-*,  ±SD).  Reduced  uptake  by  cells  on 
the  filter  was  observed  in  other  experiments. 

To  test  our  assumption  that  enhanced  activity  required  unsaturated  uptake 
kinetics,  we  repeated  the  above  flow  experiment,  but  added  radiolabeled  leucine  at  a 
concentration  expected  to  saturate  uptake  kinetics.  The  concentration  of  leucine  in  aged 
samples  for  this  experiment  was  1.1  nM,  and  was  below  detection  (<1  nM)  in 
bacteria-free  samples.  When  10.8  nM  [^H]leucine  was  added  to  the  bacteria-free  sample 
used  in  the  flow  experiments,  no  increase  in  uptake  was  observed  with  fluid  motion. 
This  confirmed  our  hypothesis  that  unsaturated  kinetics  was  a  necessary  prerequisite  for 
enhanced  microbial  uptake  with  fluid  motion. 

We  also  examined  the  effect  of  fluid  motion  on  [^HJglucose  uptake  using  aged 
samples  from  the  Roosevelt  Inlet.  We  were  unable  to  detect  any  increase  in  fluid  motion 
on  uptake.  [^H]glucose  uptake  by  suspended  averaged  3.9±0.3  pg  ml-*  min-*,  when 
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glucose  was  added  at  a  concentration  of  1  nM.  At  fluid  velocities  between  8-48  m  d'*, 
glucose  uptake  was  an  average  of  3.1±0.4  pg  ml"*  min-*,  or  about  23%  lower  than 
suspended  samples. 

At  shear  rates  of  0-2.1  s-*,  fluid  shear  did  not  appreciably  affect  leucine  uptake 
of  either  fresh  or  aged  samples  by  more  than  ±10%.  These  results  are  consistent  with 
previous  pure  culture  experiments  and  mass  transfer  calculations  for  small  molecular 
weight  compounds  with  large  diffusion  coefficients.  The  lack  of  an  effect  of  fluid  shear 
arises  from  the  inability  of  shear  to  significantly  alter  transport  rates  of  chemicals  to  the 
cell  surface  compared  to  transport  due  to  diffusion. 

The  area  covered  by  outer  membrane  proteins  involved  in  leucine  transport 
(leucine  porters)  is  important  for  examining  the  impact  of  flow  on  leucine  utilization  and 
for  understanding  how  these  bacteria  survive  in  the  low  organic  concentrations  of 
seawater.  We  estimated  that  as  much  as  9%  of  the  surface  area  of  free-living  marine 
bacteria  was  covered  by  leucine  porters.  This  percent  coverage  is  close  to  the  surface 
area  of  porins,  which  are  abundant  nonspecific  pores  or  channels  found  in 
Gram-negative  bacteila  (e.g.  E.  coll).  Assuming  that  each  E.  coli  cell  contains 
approximately  10^  porins  with  a  surface  area  of  0.9  /am^  (Nikaido  and  Vaara  1987)  and 
a  total  cell  surface  area  of  6.7  ^m^,  approximately  13%  of  the  cell  surface  is  occupied 
by  porins.  The  surface  area  of  specific  leucine  porters  in  E.  coli  is  probably  much  less 
than  for  marine  bacteria,  although  we  are  unaware  of  any  published  estimates. 

Conclusions:  Fluid  motion  can  substantially  increase  uptake  by  microorganisms, 
but  this  effect  is  a  function  of  several  factors,  including  the  substrate  type  and 
concentration.  Marine  bacteria  held  on  filters  in  a  flow  stream  of  10  m  d"^  assimilated 
leucine  8  times  faster  than  bacteria  under  similar  conditions  in  the  absence  of  fluid 
motion.  This  enhancement  was  only  observed  when  leucine  uptake  was  not  saturated. 
In  order  for  fluid  motion  to  increase  uptake,  concentration  gradients  must  be  present  at 
the  cell  surface.  At  high  leucine  concentrations  (-11  nM),  fluid  motion  had  no  effect  on 
leucine  uptake  probably  because  of  the  lack  of  a  concentration  boundary  layer.  These 
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results  imply  that  under  nutrient  limited  conditions,  attachment  to  detritus  could  increase 
bacterial  uptake  if  there  is  sufficient  advective  flow  past  the  cell.  A  detailed  description 
of  this  research  is  contained  in  Logan  and  Kirchman  (1990). 

THE  POTENTIAL  FOR  INCREASED  UPTAKE  BY  FLOCCULATING  DIATOMS 

Blooms  of  chain-forming  diatoms  commonly  flocculate  into  centimeter-sized 
aggregates  of  living,  vegetative  cells  following  nutrient  depletion  in  surface  waters  off 
southern  California.  We  examined  the  hypothesis  that  diatom  cells  within  aggregates 
experience  increased  nutrient  uptake  relative  to  unattached  cells.  We  used  data  on 
aggregates  collected  from  surface  waters  in  the  Santa  Barbara  Channel.  Three  species  of 
diatoms,  Chaetoceros  radicans.  Ch.  debilis.  and  Nitzschia  sp..  comprised  86±4%  of  the 
cells  examined.  The  remaining  cells  were  classified  as  centrics,  pennates,  or 
dinoflagellates.  Aggregates  were  very  loosely  held  together  by  tangled  spines. 

Measured  settling  velocities  of  diatom  floes  in  situ  ranged  between  49  and  190  m 
d~*,  with  an  average  of  110±40  m  d-*.  Predicted  settling  velocities  of  the  12  diatom 
floes  ranged  between  73  and  120  m  d"^  and  were  not  statistically  different  (t-test)  from 
the  observed  settling  velocities  of  these  floes,  even  though  the  calculated  velocities  did 
not  appreciably  increase  with  aggregate  diameter.  The  average  predicted  settling  velocity 
of  100±12  m  d-*  did  not  differ  significantly  from  the  in  situ  settling  velocity  of  110  in 
d'l.  Similar  values  for  the  observed  and  predicted  settling  velocities  supported  the 
calculation  of  an  equivalent  aggregate  radius  used  in  this  analysis,  and  the  application  of 
the  permeable  aggregate  model  for  diatom  floes. 

Measured  sinking  velocities  were  used  to  calculate  intra-aggregate  flow  velocities 
between  24  and  160  s-*  (2-14  m  d'*),  with  an  average  intra-aggregate  velocity  of 

70  n m  s_l  (6  m  d-^),  or  about  5%  of  the  observed  aggregate  settling  velocity  for  12 
diatom  Toes.  A  wide  range  of  intra-aggregate  velocities  of  60-160  /im  s-*  are  predicted 
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for  aggregates  in  a  narrow  size  range  of  17-20  mm  in  diameter.  This  reflects  the 
variable  porosity,  and  therefore,  permeability,  of  diatom  floes. 

For  12  diatom  floes,  we  calculated  Relative  Uptake  Factors  between  1. 0-2.1 
(average  of  1.4),  indicating  that  diatoms  within  the  aggregate  could  potentially  utilize 
dissolved  nutrients  up  to  2.1  times  faster  than  unattached  diatoms  experiencing  fluid 
shear  in  the  water  column.  This  analysis  accounts  for  nutrient  depletion  within  the 
aggregate,  and  requires  that  both  unattached  diatoms  and  the  diatom  floe  are  exposed  to 
low  nutrient  concentrations  where  a  first  order  (unsaturated)  uptake  model  accurately 
describes  nutrient  uptake  kinetics.  For  these  calculations,  we  assumed  unattached 
diatoms  existed  in  laminar  fluid  shear  (G=l  s-*). 

Relative  Uptake  Factors  were  also  calculated  assuming  unattached  cells  sank  in 
undisturbed  fluid.  Assuming  a  sinking  velocity  of  1  m  d"1  for  unattached  cells,  the 
Relative  Uptake  Factors  of  cells  within  aggregates  ranged  between  1.1  and  2.2.  This 
suggests  that  consideration  of  unattached  cells  as  either  sinking  or  suspended  within  a 
laminar  shear  field  is  unimportant  in  obtaining  Relative  Uptake  Factors  greater  than 
unity  for  diatom  floes. 

Our  data  suggest  that  the  flocculation  of  diatom  blooms  may  be  immediately 
adaptive  because  aggregation  alters  the  fluid  environment  of  a  cell.  A  single  unattached 
cell  or  chain  is  contained  within  a  microscale  eddy  and  must  move  with  the  bulk  fluid; 
however,  fluid  flow  around  and  through  a  rapidly  sinking  aggregate  actually  alters  the 
fluid  environment  of  associated  cells,  and  increases  nutrient  uptake  by  attached  cells 
compared  to  unattached  cells.  Flocculation  may,  therefore,  be  advantageous  to 
nutrient-stressed  diatoms  by  increasing  the  potential  for  uptake  of  scarce  nutrients 
following  an  intense  phytoplankton  bloom. 

Conclusions:  Using  permeability-porosity  relationships,  we  calculated 
intra-aggregate  flow  velocities  of  20-160  ixm  s'*  through  diatom  aggregates.  Although 
subject  to  considerable  uncertainty,  a  Relative  Uptake  Factor  analysis,  based  on  mass 
transfer  equations,  indicated  that  diatoms  fixed  within  aggregates  undergoing 
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gravitational  settling  could  take  up  nutrients  up  to  2.1±0.4  times  faster  than  unattached 
diatoms  experiencing  laminar  shear.  Increased  nutrient  uptake  by  aggregated  diatoms 
may  be  important  in  understanding  the  reasons  for  diatom  floe  formation.  A  detailed 
description  of  this  research  is  contained  in  Logan  and  Alldredge  (1989). 

ENHANCED  UTILIZATION  OF  '1ACROMOLECULES  BY  SUSPENDED  BACTERIA 
IN  SHEARED  FLUIDS 

The  purpose  of  this  final  aspect  of  research  was  to  investigate  the  effect  of  fluid 
shear  on  microbial  uptake  of  high  molecular  weight  dissolved  organic  compounds.  We 
examined  the  effect  of  shear  on  uptake  of  two  model  macromolecules,  an  albumin 
(65,000  amu)  as  an  example  of  a  well  characterized  protein,  and  a  dextran  (70,000  amu) 
as  an  example  of  a  complex  carbohydrate. 

Pure  cultures  of  microorganisms  were  grown  in  suspended  culture  using  a 
macromolecule  as  a  carbon  source  in  a  defined  mineral  salts  media.  We  grew  Z* 
rami&era  cultures  only  on  a  single  protein  (bovine  serum  albumin).  Since  this  culture 
was  unable  to  grow  on  dextran,  we  used  E.  coii  in  dextran  degradation  studies.  In  order 
to  prepare  cultures  for  uptake  studies,  we  processed  each  culture  through  a  series  of  3 
centrifugation/concentration  steps  (5000  g)  rinsing  the  cells  with  mineral  salts  media. 
Prior  to  uptake  experiments,  all  cell  suspensions  were  filtered  through  a  5  /im  filter  to 
remove  floes. 

We  performed  two  types  of  uptake  studies.  In  the  first  set  of  experiments,  we 
used  radiolabeled  albumin  and  dextran.  The  albumin  was  labeled  with  -^H  using 
procedures  adapted  from  Jentoft  and  Dearborn  (1979).  ^H-dextran  was  purchased  from 
SIGMA.  The  uptake  study  was  carried  out  in  60-ml  BOD  bottles,  and  uptake  was 
measured  as  filterable  radioactivity  using  0.2-jrm  polycarbonate  filters.  We  monitored 
label  incorporated  in  the  cells  in  cultures  stirred  with  magnetic  stir  bars  versus 
undisturbed  cultures.  We  found  that  uptake  of  the  label  was  2-4  times  as  large  in 
stirred  cultures  as  in  undisturbed  cultures.  In  protein  degradation  studies,  we  monitored 
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the  accumulation  of  label  in  a  <10,000  size  fraction  using  1.5  ml  centricon  centrifuge 
membrane  filters  (Millipore  Co.).  We  found  that  label  accumulated  in  the  <10,000  amu 
fraction  faster  than  in  the  undisturbed  sample,  indicating  increased  breakdown  of  the 
macromolecules  with  fluid  shear. 

In  the  second  series  of  experiments,  we  added  higher  concentrations  of  proteins 
and  dextrans  (~10  mg/1)  to  suspended  cultures,  and  monitored  dissolved  oxygen 
concentration  in  mixed  and  undisturbed  cultures.  We  found  that  the  oxygen  utilization 
rate  of  mixed  cultures  was  larger  than  undisturbed  cultures,  in  agreement  with  radiolabel 
experiments. 

An  alternative  explanation  for  the  increased  uptake  of  oxygen  and 
macromolecules  is  that  through  fluid  shear,  we  injured  cells  and  increased  cell 
respiration  rates.  However,  when  glucose  was  added  to  cultures,  uptake  of  ^H-glucose 
and  oxygen  utilization  rates  of  mixed  and  undisturbed  cells  was  not  significantly 
different.  In  addition,  the  uptake  of  ^H-leucine  was  not  affected  by  mixing,  in 
agreement  with  the  biohydrodynamical  model  that  shear  does  not  significantly  affect 
uptake  of  small  molecules  with  high  diffusion  coefficients. 

Conclusions:  Fluid  shear  can  significantly  increase  the  rate  of  macromolecule 
degradation,  as  indicated  through  enhanced  accumulation  radiolabel  in  the  cell,  enhanced 
accumulation  of  macromolecule  breakdown  products  (<10,000  amu),  and  increased 
oxygen  utilization  rates  of  mixed  cultures  versus  undisturbed  cultures.  A  detailed 
description  of  this  research  can  be  found  elsewhere  (Confer,  1990).  A  manuscript  for 
publication  in  Applied  and  Environmental  Microbiology  is  in  progress. 
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The  effect  of  fluid  flow  and  laminar  shear  on  bacterial 
uptake  was  examined  under  conditions  representative 
of  the  fluid  environment  of  unattached  and  attached 
cells  in  wastewater  treatment  bioreactors.  Laminar  shear 
rates  below  50  s'1  did  not  increase  leucine  uptake  by  sus¬ 
pended  cultures  of  Zoogloea  ramigera.  However,  leucine 
uptake  by  cells  fixed  in  a  flow  field  of  ~  1  mm  s'1  was 
55-65%  greater  than  uptake  by  suspended  cells.  En¬ 
hanced  microbial  uptake  with  advective  motion  is 
consistent  with  mass  transfer  rates  calculated  using 
Sherwood  number  correlations.  Advective  flow  increases 
microbial  uptake  by  increasing  collisions  between  sub¬ 
strate  molecules  and  cells  through  compression  of  the 
concentration  boundary  layer  surrounding  a  cell.  The 
rate  of  leucine  uptake  suggests  that  binding  proteins 
used  to  transport  leucine  into  the  cell  can  occupy  ap¬ 
proximately  1%  of  the  cell  surface  area. 


INTRODUCTION 

Microbial  adhesion  to  surfaces,  such  as  rocks  or  plastic 
media  in  trickling  filters,  and  to  suspended  inorganic  parti¬ 
cles  and  other  microorganisms  in  activated  sludge  reactors, 
is  crucial  for  the  performance  of  all  biological  wastewater 
treatment  processes.1 2  Although  the  environmental  factors 
that  induce  and  control  microbial  adhesion  are  not  well  un¬ 
derstood,  it  is  evident  that  permanent  bacterial  attachment 
requires  the  production  of  extracellular  polymers.3  4  In  a 
competitive  environment,  carbon  and  energy  requirements 
for  polymer  synthesis  must  be  less  than  the  benefits  re¬ 
ceived  from  attachment.  This  suggests  that  attachment  re¬ 
sults  in  increased  nutrient  availability. 

In  most  wastewater  treatment  bioreactors,  nutrient  con¬ 
centrations  are  low,  and  microorganisms  grow  at  rates 
substantially  less  than  their  maximum.  At  these  concentra¬ 
tions,  diffusion  of  nutrients  to  suspended  cells  may  limit 
nutrient  uptake.  Microbial  attachment  alters  the  fluid  envi¬ 
ronment  of  a  microorganism;  while  unattached  bacteria 
must  move  with  the  bulk  fluid,  attached  bacteria  can  have 
advective  flow  oast  their  surface.  Therefore,  attachment 
may  increase  nutrient  uptake  by  reducing  the  limitation  of 
transport  set  by  diffusion. 

Unattached  bacteria  in  both  engineered  and  natural  envi¬ 
ronments  are  substantially  smaller  than  the  Kolmogorov 
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microscale  of  turbulence  and  exist  within  an  environment 
characterized  as  laminar  shear.  In  bioreactors,  the  mean 
shear  rate  is  calculated  from  the  power  input  by  fluid  mix¬ 
ing  and  aeration  and  is  estimated5  to  be  in  the  range  of  90- 
220  s  ',  although  shear  rates  may  be  70  times  higher  near 
impellors.6  Natural  systems,  such  as  lakes  and  oceans,  are 
characterized  by  much  lower  shear  rates  in  the  range  of 
<  0. 1-10  s'1.1  Shear  rates  would  have  to  approach  106  s'1 
to  reduce  the  microscale  of  turbulence  to  the  size  of  a 
l-/zm  cell. 

Attached  microorganisms  exist  in  a  variety  of  environ¬ 
ments.  Microbes  can  form  aggregates  as  well  as  attach  to 
suspended  particles,  rising  bubbles,  and  fixed  surfaces. 
Biological  aggregates  are  highly  porous,  and  theoretical 
calculations8' 10  and  experimental  evidence"'13  support  fluid 
motion  through  these  aggregates  at  velocities  approaching 
the  free  field  velocity.  Floes  from  an  activated  sludge  pro¬ 
cess  0.5- 1.5  mm  in  longest  length  sink  at  rates  of  0.2-2 
mm  s'1  (ref.  14)  and  bubbles  formed  during  bioreactor 
aeration  rise  at  velocities  up  to  200-300  mm  s'1  (ref.  15) 
and  may  transport  both  cells  and  floes.  Marine  snow  parti¬ 
cles,  large  amorphous  aggregates  that  form  in  marine  sys¬ 
tems  comprised  of  phytoplankton,  bacteria,  and  other 
suspended  organic  matter,  sink  at  velocities  up  to  4.5  mm 
s'1. 16  Fluid  motion  has  been  shown  to  increase  nutrient  up¬ 
take  of  phytoplankton,17'19  and  Logan  and  Hunt8  9  have  cal¬ 
culated  that  substantial  fluid  flow  through  highly  porous 
aggregates  can  result  in  an  overall  increase  in  uptake  kinet¬ 
ics  of  attached  bacteria. 

The  purpose  of  this  investigation  is  to  quantify  the  effect 
of  fluid  flow  and  fluid  shear  on  bacterial  kinetics  since 
these  two  fluid  environments  represent  the  range  of  fluid 
motions  experienced  by  attached  and  unattached  bacteria. 
In  this  study,  the  effect  of  fluid  environment  on  microbial 
growth  was  calculated  by  comparing  leucine  uptake  by 
bacteria  suspended  in  undisturbed  fluid  with  uptake  by 
bacteria  fixed  in  a  uniform  flow  field  and  suspended  in  a 
laminar  shear  field.  Leucine  uptake  was  monitored  since 
the  rate  of  leucine  incorporation  by  bacteria  has  been  de¬ 
termined  to  be  an  index  of  protein  synthesis  for  laboratory 
and  natural  assemblages  of  bacteria.20  The  measured  effect 
of  fluid  motion  on  bacterial  kinetics  was  compared  to  a 
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model  based  on  an  analysis  of  mass  transfer  to  spheres  in 
advective  flow  and  laminar  shear  fields. 

EXPERIMENTAL  METHODS 
Culture  Conditions 

Pure  cultures  of  Zoogloea  ramigera  (ATCC  1100),  a  floc- 
culant  gram  negative  rod  that  attaches  through  production  of 
cellulose  fibrils,  were  used  in  all  experiments.  Z.  ramigera 
is  thought  to  be  important  in  the  formation  of  rapidly  set¬ 
tling  floes  in  activated  sludge  wastewater  treatment  reac¬ 
tors. :i  Long-term  cultures  were  maintained  on  nutrient 
broth  agar  (Difco)  and  transferred  to  sterile  liquid  media- 
for  each  experiment.  Liquid  media  consisted  of  1  g  L"1  of 
glucose  in  a  mineral  salt  buffer  (MSB),  containing,  per  liter 
of  water:  0.57  g  NH4NO„  0.5  g  KH2P04,  1  g  K,HP04, 
0.2  g  MgS04,  0.04  mg  FeCl3,  and  0.02  mg  vitamin  B- 12. 
All  cells  used  in  these  experiments  were  grown  in  a  shaker 
table,  harvested  during  exponential  growth,  and  prefiltered 
(to  remove  floes)  through  a  8-/xm  cellulose  acetate  filter 
(Millipore  Corp.)  at  100  mm  Hg  immediately  before  use. 
All  studies  were  performed  within  several  hours  using 
samples  isolated  from  the  same  suspended  culture. 

Leucine  uptake  by  cells  of  Z.  ramigera  was  determined 
using  radiolabeled  leucine  (4,5-[JH]leucine,  60  Ci  mmol"1, 
ICN  Laboratories),  supplemented  with  unlabeled  leucine, 
at  final  concentrations  of  10-3— 102  mg  L'1.  All  samples 
were  analyzed  using  a  Beckman  5801  liquid  scintillation 
counter  with  H-number  correction  for  sample  quenching. 
Microorganisms  were  enumerated  using  acridine  orange- 
epifluorescence  direct  count  methods.22  Cells  were  ap¬ 
proximately  0.8  fi.m  wide  by  2  /im  long,  equivalent  in 
volume  to  a  spherical  cell  with  0.7  /im  radius. 

Effect  of  Fluid  Flow 

The  effect  of  fluid  motion  on  bacterial  uptake  was  mea¬ 
sured  by  comparing  leucine  uptake  of  suspended  cells  to 
leucine  uptake  by  cells  fixed  in  a  flow  field.  Cells  were 
maintained  in  a  flow  field  by  vacuum  filtration  onto  a 
25-mm  filter  (various  types),  supported  by  a  5-jzm  cellu¬ 
lose  acetate  backing  filter.  These  filters  were  supported  by 
stainless  steel  bases  and  funnels  in  a  10-place  vacuum  box 
fitted  with  Teflon  valves  (Hoefer  Scientific).  Three  differ¬ 
ent  types  of  filters  were  used  to  capture  bacteria:  0.2  gm 
polycarbonate  filters  (Poretics  Corp.),  0.2  ptm  aluminum 
oxide  filters  (Anotech  Corp.),  and  0.45  (im  cellulose  ace¬ 
tate  filters  (Gelman  Corp.).  The  0.2-/zm  filters  are  surface 
filters,  and  cells  were  completely  retained  on  the  surface  of 
the  filter.  The  0.45-jum  filter  is  a  depth  filter,  and  some 
penetration  of  cells  into  the  filter  matrix  may  have  occurred. 

Cells  were  fixed  in  a  flow  field  by  vacuum  filtration  of  a 
solution  consisting  of  200- /xL  of  prefiltered  cells  with  an 
absorbance  of  0.42  in  2  mL  of  MSB.  The  filter  was  imme¬ 
diately  rinsed  with  2  mL  of  sterile  MSB  until  only  a  light 
film  of  MSB  solution  remained  on  the  filter  (to  prevent 
cell  desiccation).  Ceil  loadings  on  the  filter,  by  projected 


surface  area,  were  less  than  1%.  A  solution  of  sterile  MSB 
containing  ~  0.5  /xg  L"1  [’HJleucine  was  pulled  through 
the  filter  at  different  velocities  by  varying  the  vacuum.  In¬ 
cubation  times  between  1  and  2  min  were  achieved  by 
varying  flow  volumes.  After  fluid  had  been  pulled  through 
the  filter,  bacteria  were  immediately  inactivated  using 
2  mL  of  2%  formalin.  Filters  were  rinsed  with  an  addi¬ 
tional  5  mL  of  sterile  MSB,  withdrawn  from  the  vacuum 
box,  and  combined  with  10  mL  of  ACS  scintillation  cock¬ 
tail  (Amersham  Corp.)  for  counting.  In  separate  experi¬ 
ments,  repeated  washing  of  cells  with  fresh  media  did  not 
reduce  radiolabel  concentration,  indicating  formalin  addi¬ 
tion  did  not  lyse  cells.  Further  addition  of  radiolabeled 
leucine  after  formalin  addition  resulted  in  no  leucine  up¬ 
take,  indicating  a  loss  of  cell  activity. 

Uptake  by  suspended  cells  was  determined  by  incuba¬ 
tion  of  200  fiL  of  prefiltered  cells  with  2  mL  of  MSB  con¬ 
taining  [3H]leucine  for  same  length  of  time  as  the  flow 
experiment.  After  incubation,  cells  were  inactivated  using 
formalin,  and  the  solution  was  filtered,  rinsed,  and  ana¬ 
lyzed  as  described  above.  Both  the  flow  and  suspended- 
counts  were  corrected  by  subtraction  of  abiotic  blanks 
prepared  under  identical  conditions. 


Effect  of  Fluid  Shear 


A  defined  laminar  shear  environment  was  obtained  using  a 
rotating  cylinder  apparatus.23  The  device  consisted  of  two 
vertically  aligned  concentric  acrylic  cylinders  separated  by 
a  9-mm  gap.  The  fixed  inner  cylinder  was  8  cm  in  diame¬ 
ter  and  15  cm  high.  The  outer  cylinder  rotated  at  a  fixed 
speed  producing  a  uniform  shear  rate  between  the  cylin¬ 
ders.  The  shear  rate  G  was  calculated23  using 


2its  2R,R2 
"60  (*!-*?) 


(1) 


where  s  is  the  rotational  speed  (rpm)  and  /?,  and  R2  are  the 
inner  and  outer  radii,  respectively. 

Prefiltered  cells  (2  mL)  were  combined  with  500  mL  of 
MSB  containing  [3H]Ieucine,  and  400  mL  of  this  solution 
was  poured  into  the  rotating  cylinder  annulus.  The  cylinder 
was  started,  and  time  zero  samples  were  obtained  by  with¬ 
drawing  three  10-mL  samples  from  the  remaining  undis¬ 
turbed  100-mL  suspended-cell  solution.  After  15  min,  three 
10-mL  samples  were  withdrawn  from  the  rotating  cylin¬ 
der  solution  at  a  depth  of  5  cm  and  from  the  undisturbed 
suspended-cell  solution.  Samples  were  combined  with  for¬ 
malin  solution  (2%  final  concentration)  and  filtered  using 
0.2-^im  polycarbonate  filters.  Uptake  was  calculated  as  the 
difference  between  initial  and  final  radioactivity  of  the  filters. 


KINETIC  MODEL 


Assuming  Monod  kinetics,  the  change  in  nutrient  mass 
concentration  with  time,  dC/dt,  for  a  microbial  culture  is 


HrfC  _  jiC, 

N  dt  Y{K  +  C.) 


(2) 
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where  N  is  the  cell  number  concentration,  C  and  C ,  are  the 
bulk  and  cell  surface  nutrient  concentrations,  respectively, 
jx  is  the  maximum  growth  constant,  k  is  the  half  saturation 
constant,  and  Y  is  a  yield  coefficient.  Since  Cs  is  not  easily 
determined,  the  kinetic  constants  are  usually  calculated  as¬ 
suming  C  =  C,.24  At  low  nutrient  concentrations,  C  K, 
and  Equation  (2)  can  be  simplified  to  a  first-order  kinetic 
expression. 


^IdC  _  £C 
N  dt  ~  YK 


(3) 


As  written.  Equation  (3)  is  independent  of  fluid  motion. 
The  effect  of  fluid  motion  on  microbial  kinetics  can  be  de¬ 
termined  by  describing  nutrient  utilization  using  a  mass 
transfer  analysis. sv  For  a  suspended  culture  containing  only 
free-living  cells,  the  rate  of  nutrient  utilization  shown  in 
Equation  (2)  is  equal  to  the  rate  of  mass  transfer  to  each 
cell,  Q,  or 

=  Q  =  kW(C  -  C,)  (4) 

N  dt 


where  k  is  a  mass  transfer  coefficient  and  a  is  the  cell 
radius.  Fluid  motion  increases  microbial  uptake  through 
compression  of  the  concentration  boundary  layer  surround¬ 
ing  the  cell. 

The  mass  transfer  coefficient  can  be  written  in  terms  of 
the  dimensionless  Sherwood  number,  Sh  =  k a/D,  where 
D  is  the  nutrient  diffusivity.  For  a  spherical  cell  in  stagnant 
fluid,  the  Sherwood  number  is  unity.  Incorporating  the 
Sherwood  number  into  Equation  (3)  and  assuming  C,  is 
very  small, 

Q  =  4t toD  Sh  C  (5) 

Several  Sherwood  number  correlations  have  been  pro¬ 
posed  to  account  for  the  effect  of  fluid  environment  on  the 
rate  of  mass  transfer  to  spheres.  At  Reynolds  numbers 
much  less  than  unity  (Re  =  ua/v  <  1),  Brian  and  Hales 25 
correlated  the  Sherwood  number,  Sh„ ,  for  a  sphere  fixed 
in  a  flow  stream  of  velocity  u,  with  the  Peclet  number, 
Pe„  =  ua/D  <  104,  as 

Shu  “  (1  +  0.48  Pef)1'2  (6) 

Munk  and  Riley26  determined  the  effect  of  sinking  velocity 
on  mass  transfer  to  phytoplankton  using 

Sh„  =  1  +  0.5  Pe  +  0.6  Pe2  (7) 

The  effect  of  fluid  velocity  on  phosphorous  uptake  by  the 
diatom  Thalassiosira  fluviatilis,  determined  by  Canelli  and 
Fuhs,18  was  used  by  Logan  and  Alldredge10  to  calculate  the 
Sherwood  number  in  water  as 

Sh,  =  1.45  Re0  57  (8) 

A  Sherwood  number  for  a  sphere  in  laminar  shear 
flow,  Sh0 ,  was  derived  by  Frankel  and  Acrivos  for  Pe  = 
axG /D  «  1  as 

Shc  =  1  +  0.26  Pe05  (9) 


Since  the  Sherwood  number  is  unity  for  a  cell  in  stagnant 
fluid,  the  ratio  of  uptake  by  a  cell  in  a  fluid  field  or  in  lami¬ 
nar  shear,  Q,  to  uptake  by  a  cell  in  undisturbed  or  stagnant 
fluid,  Q.,  calculated  using  Equation  (5),  is 

^-  =  Sh  (10) 

Therefore,  the  ratio  of  mass  transfer  to  cells  in  these  differ¬ 
ing  fluid  environments  to  cells  in  stagnant  fluid  conditions 
can  be  determined  using  one  of  the  Sherwood  number  cor¬ 
relations  in  Equations  (6-9). 

RESULTS 

Suspended  Cultures 

The  kinetic  models  presented  in  Equations  (3-6)  require 
nutrient  uptake  by  microorganisms  to  be  a  linear  function 
of  concentration.  We  examined  leucine  uptake  by  suspended 
cultures  of  Z.  ramigera  at  three  different  incubation  times, 
at  concentrations  between  0.5  fig  L“‘  and  100  mg  L"1. 
As  shown  in  Figure  1 ,  leucine  uptake  was  linear  with  con¬ 
centration  over  a  range  spanning  six  orders  of  magnitude. 
As  incubation  times  increased,  the  rate  of  uptake  decreased, 
suggesting  feedback  from  an  intracellular  leucine  pool.  For 
incubation  times  of  0.5,  1 ,  and  7  min,  the  slopes  ( ±  standard 
error)  of  the  lines  shown  in  Figure  1  were  7.0  ±  0.02, 
6.4  ±  0.02,  and  5.2  ±  0.03  x  10~3  min'1,  respectively. 
The  r2  values  for  all  lines  were  greater  than  0.98.  These  re¬ 
sults  indicate  that  saturation  kinetics  were  not  exhibited  by 
Z.  ramigera  over  a  wide  range  of  leucine  concentrations. 
A  similar  lack  of  saturation  kinetics  over  a  wide  range  of 
substrate  concentrations  was  also  observed  by  Rubin  and 
Alexander2*  for  phenol  mineralization  by  natural  bacterial 
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LOG  LEUCINE  (mg/I) 

Figure  1.  Leucine  uptake  as  a  function  of  leucine  concentration,  at  differ¬ 
ent  incubation  times:  (□)  0.5  min,  (x)  1.0  min.  and  (A)  7  min.  All  lines 
have  an  rz  value  greater  than  0.98. 
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assemblages  at  phenol  concentrations  between  0.002  and 
200  jug  L” 1 . 

Effect  of  Fluid  Flow 

Several  flow  experiments  conducted  with  Z.  ramigera  held 
on  0.2-jum  polycarbonate  filters  indicated  leucine  uptake  by 
cells  fixed  in  a  flow  field  was  55-65%  greater  than  leucine 
uptake  by  suspended  cells.  Shown  in  Figure  2  are  two  rep¬ 
resentative  fluid  flow  experimental  results.  In  the  first  ex¬ 
periment  (Fig.  2a),  only  radiolabeled  leucine  was  added  to 
the  culture.  In  the  second  experiment  (Fig.  2b)  cells  were 
acclimated  to  10  ng  L'1  of  leucine  for  30  min  prior  to  addi¬ 
tion  of  the  radiolabel  and  the  start  of  the  experiment.  In 


these  experiments,  leucine  uptake  increased  with  fluid  ve¬ 
locity  at  velocities  less  than  ~1  mm  s’1.  Little  or  no  in¬ 
crease  in  leucine  uptake  was  observed  at  higher  velocities 
between  1  and  2.5  mm  S'1.  On  occasion,  some  cells  were 
completely  inactivated  at  higher  flow  velocities  and  did  not 
take  up  any  leucine,  probably  due  to  cell  disiuption  at 
vacuums  of  600  mm  Hg  needed  to  obtain  higher  velocities. 
Fluid  velocities  below  1  mm  s’1  were  obtained  at  vacuums 
below  250  mm  Hg,  a  range  indicated  by  Fuhrman  and 
Bell29  to  prevent  significant  release  of  amino  acids  from 
marine  microorganisms. 

Fluid  flow  experiments  using  0.45  cellulose  acetate  and 
0.2-fj.m  aluminum  oxide  filters  to  hold  bacteria  were  char¬ 
acterized  by  large  variations  in  bacterial  uptake  about 


0.0  0.5  1.0  1.5  2.0  2.5 

VELOCITY  (mm/s) 

(a) 


Figure  2.  The  effect  of  fluid  velocity  on  leucine  uptake  for  cells  held 
on  a  0.2-^im  polycarbonate  filter  at  different  leucine  concentrations: 
(a)  0.4  fig  L"1;  (b)  10.8  fig  L~'.  Lines  are  fitted  by  eye.  Error  bars  corre¬ 
spond  to  ±  l  SD  of  3-5  measurements. 
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mean  values  (Fig.  3).  Cellulose  acetate  filters  retained 
large  concentrations  of  ['Hjleuei.ic  label  during  flow  ex¬ 
periments  which  could  not  be  removed  by  filter  rinsing. 
This  substantially  increased  the  amount  of  abiotic  adsorp¬ 
tion  of  leucine  to  filters.  These  increased  background 
leucine  concentrations  produced  large  variations  in  uptake 
for  a  given  fluid  velocity.  Based  on  averages  shown  in  Fig¬ 
ure  3a,  the  effect  of  fluid  motion  on  uptake  was  similar  to 
that  observed  with  the  polycarbonate  (surface)  filters.  A 
maximum  average  increase  of  60^  was  consistent  with  pre¬ 
vious  results  using  the  polycarbonate  filters.  However,  the 
larger  errors  associated  each  measurement  did  not  result  in 
a  significant  correlation  of  fluid  velocity  with  uptake. 


The  results  shown  in  Figure  3b  for  the  aluminum  oxide 
Filter  suggest  an  increase  in  uptake  with  fluid  flow,  but  the 
large  variability  of  the  data  also  do  not  present  any  signifi¬ 
cant  pattern.  Flow  experiments  with  the  aluminum  oxide 
filters  required  much  higher  vacuums  (>250  mm  Hg)  to 
achieve  the  same  fluid  velocities  used  in  polycarbonate  fil¬ 
ter  experiments.  These  high  vacuums  likely  impaired  net 
amino  acid  uptake. 

The  effect  of  fluid  motion,  predicted  from  a  mass  trans¬ 
fer  analysis,  is  shown  in  Figure  4.  Sherwood  numbers 
were  calculated  using  Equations  (6-8),  an  average  cell 
size  of  0.7  /am,  and  a  molecular  diffusivity  for  leucine  of 
0.7  x  10" 5  cm:  s  '  determined  using  the  Wilke-Chang 


Figure  3.  The  effect  of  fluid  velocity  on  leucine  uptake  foe  cells  held  on: 
(a)  0.43  /am  cellulose  acetate  filter  at  0.3  /ig  L1  of  leucine;  (b)  0.2  4m  alu¬ 
minum  oxide  filter  at  0.7  fig  L"1  of  leucine.  Error  bars  correspond  to  ±1  SD 
of  3-5  measurements 
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Figure  4.  The  effect  of  fluid  velocity  presented  in  terms  of  a  Sherwood 
number  Data  from  Fig.  2  are  plotted  as  a  ratio  of  uptake  by  cells  in  a  flow 
field  to  cells  in  stagnant  fluid.  Lines  are  plotted  using  Sherwood  number  cor¬ 
relations  in  Eqs.  (6)-(8). 


correlation.’0  These  calculations  were  based  on  the  as¬ 
sumption  that  leucine  concentrations  at  the  cell  surface 
were  much  less  than  the  bulk  concentration,  that  uptake 
was  not  zero  order,  and  that  flow  past  cells  was  completely 
uniform.  At  a  flow  velocity  of  1  mm  s’1.  Equations  (6), 
(7),  and  (8)  predict  an  increase  in  mass  transfer  of  22,  110, 
and  65%,  respectively.  It  can  be  seen  in  Figure  4  that  the 
correlation  used  by  Munk  and  Riley26  consistently  overes¬ 
timated  the  effect  of  fluid  motion.  The  mass  transfer  corre¬ 
lation  developed  by  Logan  and  Alldredge'0  [Equation  (8)] 
gave  the  best  fit  to  the  experimental  results  of  this  study. 

Effect  of  Fluid  Shear 

The  effect  of  fluid  shear  on  leucine  uptake  by  Z.  ramigera 
is  shown  in  Figure  5  as  a  ratio  of  uptake  by  sheared  cells 
to  cells  fixed  in  stagnant  fluid.  Although  results  at  fluid 
shear  rates  above  50  s'1  are  shown  in  Figure  5b,  the  fluid 
within  the  cylinder  became  unstable  at  these  shear  rates 
and  could  no  longer  be  considered  laminar.  Values  at  zero 
shear  were  determined  by  not  turning  on  the  rotating  cylin¬ 
der  after  sample  addition.  This  “zero  shear”  case  was  only 
within  20%  of  values  obtained  from  analysis  of  suspended 
cultures.  Thus,  the  experimental  error  of  shear  experiments 
is  very  large  and  at  least  20%.  Within  this  large  experi¬ 
mental  error,  fluid  shear  does  not  increase  leucine  uptake 
by  bacteria.  Moreover,  there  is  a  consistent  trend  of  de¬ 
creasing  uptake  with  shear  that  suggests  that  higher  lami¬ 
nar  shear  rates  impaired  leucine  uptake  by  suspended 
bacteria. 

A  mass  transfer  analysis  indicates  that  fluid  shear  should 
not  increase  leucine  uptake  by  more  than  5%  at  shear  rates 
below  50  s'1.  This  is  shown  in  Figure  6  using  Equations 
(9)  and  (10)  compared  to  data  from  Figure  5.  The  mathe¬ 


matical  analysis  is  limited  to  Peclet  numbers  much  less 
than  unity.  For  a  cell  radius  of  0.7  /zm  and  leucine  diffu- 
sivity  of  0.7  x  10'5  cm2  s'1,  a  Peclet  number  of  0.1  oc¬ 
curs  at  a  shear  rate  of  140  s'1  and  is  larger  than  shear  rates 
examined  in  this  study. 

DISCUSSION 

From  both  experimental  results  and  theoretical  consider¬ 
ations,  fluid  shear  is  not  expected  to  substantially  increase 
microbial  uptake  by  unattached  bacteria.  Uptake  experi¬ 
ments  with  Z.  ramigera  showed  that  fluid  shear  and,  there¬ 
fore,  fluid  mixing  intensity,  did  not  increase  leucine  uptake 
by  bacteria  suspended  in  a  laminar  shear  field.  Calcula¬ 
tions  based  on  a  Sherwood  number  correlation  [Equa¬ 
tion  (6)]  predicted  a  maximum  increase  of  only  5%  at  shear 
rates  below  50  s'1.  This  calculation  is  also  in  agreement 
with  the  prediction  of  Purcell31  that  fluid  shear  at  intensi¬ 
ties  typical  of  bioreactors  could  not  substantially  increase 
mass  transfer  of  molecules  to  bacteria-sized  particles. 

Fluid  flow  past  bacteria,  however,  can  substantially  in¬ 
crease  bacterial  uptake.  Leucine  uptake  by  Z.  ramigera 
increased  by  an  average  of  55-65%  at  flow  rates  up  to 
1  mm  s'1.  In  order  for  cells  to  be  present  in  this  type  of 
flow  field  under  natural  conditions,  they  must  either  be  at¬ 
tached  to  a  surface  or  attached  to  other  cells  in  permeable 
aggregates.  Enhanced  chemical  reaction  rates  with  fluid 
motion  have  been  demonstrated  for  forced  air  convection 
in  a  porous  catalyst  packed  bed  reactor32,33  and  for  hydro¬ 
lysis  of  benzoyl  L-arginine  ethyl  ester  by  CM-cellulose- 
fictin  enzyme  in  a  packed  bed  immobilized  enzyme 
reactor.34  Since  we  were  able  to  show  increased  bacterial 
uptake  with  fluid  motion,  our  experiments  support  predic¬ 
tions  of  enhanced  uptake  by  cells  within  permeable  aggre- 
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Figure  5.  The  ratio  of  leucine  uptake  by  cells  in  a  laminar  shear  field  to  up¬ 
take  by  cells  in  undisturbed  fluid  at  (a)  low  and  (b)  high  shear  rates.  Shear 
rates  above  50  s“ 1  are  no  longer  laminar. 


gates  for  bacterial  aggregates  suspended  in  bioreactors  and 
attached  to  bubbles.’ 

The  proposed  mass  transfer  model  of  bacterial  uptake  is 
analogous  to  a  model  proposed  by  Yao  et  al.35  describing 
the  filtration  of  molecules  by  particles  (sand  grains)  in  a 
water  filter.  The  microorganism,  like  the  inorganic  filter 
particle,  is  a  collector  of  material  transported  to  the  collec¬ 
tor  surface  by  Brownian  diffusion.  If  the  rate  of  molecular 
diffusion  limits  capture  by  the  collector,  then  the  rate  of 
capture  can  be  increased  by  advective  fluid  flow  past  the 
suiface.  This  increase  is  attributed  to  a  compression  of  the 
concentration  boundary  layer  surrounding  the  collector 
which  increases  the  mass  flux  to  the  collector  surface. 


The  filtration  model  of  Yao  et  al.”  which  uses  a  mass 
transfer  analysis  developed  by  Levich,36  is  based  on  evi¬ 
dence  that  molecules  may  collide  many  times  with  a  sur¬ 
face  before  achieving  a  successful  collision.  The  effect  of 
unsuccessful  collisions  is  incorporated  in  the  filtration 
model  using  a  parameter  defined  as  a  collector  efficiency, 
composed  of  a  probability  of  hydrodynamic  encounter  and 
a  collision  efficiency.  An  analogous  argument  is  made  here 
for  the  incorporation  of  a  collector  efficiency  parameter 
into  a  microbial  uptake/mass  transfer  model.  According  to 
the  Michaelis-Menten  description  of  enzyme  kinetics,  a 
molecule  reversibly  and  repeatedly  binds  to  an  enzyme  be¬ 
fore  successful  product  production.  Very  few  enzymes 
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Figure  6.  The  elfecl  ol  ■.hear  presented  in  terms  of  a  Sherwood  number 
Data  trom  Fie  5.  without  zero  values,  are  plotted  as  a  function  of  the  log 
shear  rate  ( ne  line  is  plotted  using  the  Sherwood  number  correlation  in 
Hq.  id i 


achieve  "kinetic  perfection."'  or  operate  at  the  limit  of 
diffusion-controlled  uptake  rates.  Therefore,  the  success 
of  repeated  molecule-enzyme  binding,  or  collision  effi¬ 
ciency.  must  be  less  than  unity  for  microbial  systems. 

A  collector  efficiency  can  be  calculated  for  leucine  up¬ 
take  using  the  experimental  results  of  this  study  for  sus¬ 
pended  cultures  of  Z.  ramigera.  During  leucine  uptake,  the 
bulk  and  surface  concentrations  of  nutrient  are  related  by 

C\  =  ( I  -  E„  )C  (II) 

where  £„,  the  collector  efficiency  of  a  cell,  is  assumed  to 
be  a  function  of  the  extent  of  molecule-enzyme  affinity  as 
well  as  the  shape  of  the  cell.  Substituting  Equation  (11) 
into  Equation  (4)  and  using  data  in  Figure  2a,  the  resulting 
collector  efficiency  is  approximately  equal  to  0.01 .  By  this 
calculation,  only  once  during  every  100  collisions  of  a 
molecule  with  the  cell  is  leucine  incorporated  into  the  cell. 

The  low  value  of  the  collector  efficiency  (EB  4  I)  im¬ 
plies  the  concentration  gradient  surrounding  a  cell  is  negli¬ 
gible  and  that  leucine  uptake  is  kinetically  controlled.  If 
this  were  true,  compression  of  the  concentration  boundary 
layer  under  these  conditions  would  not  increase  bacterial 
uptake  of  leucine  by  more  than  1%.  However,  this  as¬ 
sumes  complete  coverage  of  the  bacterial  surface  by  bind¬ 
ing  proteins  used  for  leucine  transport  across  the  bacterial 
membrane.  Complete  enzyme  coverage  on  the  cell  surface 
for  leucne  transport  is  unlikely.  The  effect  of  incomplete 
enzyme  coverage  on  mass  transfer  can  be  calculated  using 
the  analysis  of  Dudukovic  and  Mills. 18  For  their  analysis, 
the  size  of  an  enzyme  must  be  known  and  is  estimated  here 
as  0.003  /am  from  size-molecular  weight  values  given  in 
Perry  and  Chilton30  and  assuming  an  enzyme  molecular 
weight  of  100.000.  If  the  concentration  boundary  layer 
spans  0  01  /urn,  which  would  occur  for  a  1%  change  in 


concentration  across  a  l-/um  cell,  then  Dudukovic  and 
Mill's  analysis  indicates  the  mass  transfer  coefficient  must 
be  reduced  by  the  fraction  of  area  available  for  leucine  up¬ 
take.  For  example,  if  enzyme  coverage  is  only  1%  of  the 
surface  area  and  fluid  motion  reduces  the  concentration 
boundary  layer  to  within  10%  of  the  thickness  under  stag¬ 
nant  flow  conditions,  uptake  by  a  cell  fixed  in  a  flow  field 
will  be  3.3  times  as  large  as  leucine  uptake  under  stagnant 
flow  conditions.  Incomplete  enzyme  coverage  of  the  cell 
can  therefore  explain  the  observed  effect  of  fluid  motion 
on  bacterial  uptake  of  leucine. 

Other  factors  may  have  contributed  to  the  observed  ef¬ 
fect  of  fluid  motion  on  microbial  uptake,  and  -hree  addi¬ 
tional  explanations  were  examined.  First,  incomplete  label 
dissolution  during  suspended-cell  experiments  could  re¬ 
duce  leucine  uptake.  When  Rhodamine  dye  was  substi¬ 
tuted  for  the  added  radiolabel,  dissolution  of  the  dye  was 
nearly  instantaneous,  indicating  sufficient  mixing.  Second, 
the  vacuum  pressures  might  have  indirectly  increased 
leucine  uptake.  However,  very  high  vacuums  (—600  mm 
Hg)  frequently  inhibited  leucine  uptake  by  cells  and  indi¬ 
cated  a  negative  response  of  cells  to  vacuum.  Other  studies 
have  shown  high  shear  will  inhibit  the  uptake  of  another 
amino  acid  by  attached  bacteria.  Eighmy  and  Bishop39 
found  that  binding  proteins  used  for  aspartate  uptake  by 
mixed  cultuies  biofilms  were  inhibited  at  shear  rates  of 
3.1  N  m':.  Third,  there  may  have  been  a  slime  capsule 
around  the  cell  that  reduced  leucine  transport  to  the  cell 
surface.  However,  no  noticeable  capsule  was  observed  on 
cells  prefiltered  through  the  8-/xm  filter  and  used  in  these 
experiments. 

Additional  studies  reported  in  the  literature  support  the 
effect  of  enhanced  uptake  by  attached  bacteria  and  cells 
fixed  in  a  flow  field.  Kirchman  and  Mitchell40  and  Simon41 
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have  shown  that  bacteria  attached  to  particle*;  in  natural 
waters  such  as  lakes  ire  larger  and  have  higher  uptake 
rates  of  amino  acids  than  suspended  bacteria.  These  stud¬ 
ies  were  not  well  controlled  with  respect  to  tluid  hydrody¬ 
namics  and  alternative  explanations,  such  as  colonization 
of  particles  by  different  species  than  those  remaining  unat¬ 
tached.  and  increased  concentrations  of  amino  acids  ad¬ 
sorbed  to  particle  surfaces,  etc.,  could  also  explain  their 
observations.  However,  a  study  by  Canelli  and  Fuhs  *  con¬ 
clusively  demonstrated  the  effect  of  fluid  motion  on  phos¬ 
phate  uptake  by  pure  cultures  of  the  diatom  Thallasiosira 
thauirihs  held  on  a  filter  surface.  They  observed  a  130% 
increase  in  phosphorous  uptake  with  fluid  velocities  up  to 
0.4  mm  s  and  no  further  increase  in  uptake  above  this 
velocity,  based  on  their  data,  a  collector  efficiency  of 
£„  =  0.025  was  calculated  for  phosphate  uptake  by  this 
diatom.'  Thus,  an  increase  in  uptake  kinetics,  similar  to 
that  observed  in  this  study  for  uptake  of  an  amino  acid  by 
bacteria,  was  observed  by  Canelli  and  Fuhs  for  a  different 
cell  and  substrate  than  were  examined  in  this  study. 

The  results  of  the  present  study  suggest  that  the  effect  of 
tluid  motion  may  be  adequately  described  by  a  dimension¬ 
less  Sherwood  number  correlation.  The  collision  of  a 
molecule  moving  by  Brownian  motion  with  a  collector,  de¬ 
fined  in  the  filtration  model  of  Yao  et  al."  by  the  product 
of  two  terms,  the  probability  of  collision  and  the  collision 
or  sticking  efficiency,  are  combined  in  this  discussion  in 
the  constant  £„ .  the  collector  efficiency.  For  microorgan¬ 
isms.  incomplete  coverage  of  the  cell  surface  by  enzymes 
reduces  the  area  available  for  mass  transfer  to  the  cell. 
Therefore,  the  difference  in  values  of  £„  =  0.01  calcu¬ 
lated  for  leucine  uptake  by  Z.  rarniRera  and  £fl  =  0.025 
determined  for  phosphate  uptake  by  T.  Jluvuitilis  likely 
occurs  due  to  the  different  affinities  of  the  nutrient  for  the 
respective  enzyme  as  well  as  incomplete  surface  coverage 
of  enzyme.  Advective  flow  past  the  cell  either  modifies 
the  sticking  efficiency  of  the  molecule  or,  more  likely, 
enhances  the  probability  of  collisions  of  nutrient  with  the 
cell  surface. 

Based  on  the  results  of  this  study,  a  model  describing 
the  rate  of  uptake  by  a  cell  as  a  function  of  fluid  environ¬ 
ment  can  be  constructed  by  combining  Equations  (4)  and 
(9)  to  obtain 

Q  =  AnaDEgC,  Sh  (12) 

In  this  model,  the  collision  efficiency  of  the  nutrient  with 
the  cell  surface  is  limited  toO  £  £,  i  1 .  The  Sher¬ 
wood  number  has  no  physical  limits,  but  high  fluid  veloci¬ 
ties  past  the  cell  will  ultimately  damage  the  cell  and 
decrease  nutrient  uptake.  In  our  experiments,  we  typically 
observed  a  maximum  increase  of  55-65%.  although  fur¬ 
ther  increases  were  limited  by  cell  disruption  at  higher 
vacuums  needed  to  provide  higher  fluid  velocities.  Mi¬ 
croorganisms  residing  on  the  surface  of  pipes  may  reside 
in  fluid  with  much  higher  velocities  than  examined  here.12 

An  increase  in  microbial  uptake  with  fluid  motion  can 
be  an  important  competitive  advantage  for  attached  cells  in 


biological  wastewater  treatment  reactors.  In  activated 
sludge  reactors,  long  detention  times  of  cells  imply  low 
growth  rates  but  do  not  require  low  uptake  rates  of  easily 
degraded  soluble  organics.  Under  optimal  conditions,  at¬ 
tached  microorganisms  could  utilize  easily  degraded  or¬ 
ganics  up  to  65%  faster  than  unattached  cells,  assuming  all 
cells  had  the  same  enzyme  to  assimilate  the  substrate.  As  a 
result,  attached  cells  would  outcompete  unattached  cells 
for  growth  substrates.  The  effect  of  tluid  flow  on  microbial 
kinetics,  therefore,  may  serve  as  a  basis  for  understanding 
the  advantages  of  microbial  attachment  and  may  help  ex¬ 
plain  why  most  microorganisms  are  inherently  capable  of 
attachment. 
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ABSTRACT 


Fluid  motion  can  increase  the  assimilation  of  dissolved  organics  by  attached 
compared  to  free-living  microorganinsms.  This  increase  results  from  compression  of  the 
concentration  boundary  layer  surrounding  the  ceil,  and  only  occurs  when  uptake  systems 
are  not  saturated.  We  examined  the  effect  of  advective  flow  and  fluid  shear  on  the 
uptake  of  model  compounds  (leucine  and  glucose)  by  natural  assemblages  of 
heterotrophic  bacteria.  We  found  that  [^Hjleucine  uptake  by  cells  held  in  fluid  moving 
at  20-70  m  d_1  was  8  times  larger  than  uptake  by  cells  at  a  velocity  of  3  m  d~*.  This 
effect  was  only  observed  at  low  leucine  concentrations  (-1  nM),  when  uptake  was  likely 
not  saturated.  When  we  added  leucine  at  concentrations  expected  to  saturate  leucine 
uptake  (-11  nM),  fluid  motion  past  cells  did  not  affect  uptake.  Fluid  flow  past  bacteria 
did  not  increase  [^H]glucose  uptake,  and  laminar  shear  rates  of  0.5-2  s'1  did  not 
increase  either  glucose  or  leucine  uptake  by  suspended  bacteria.  These  results  indicate 
that  fluid  motion  increases  bacterial  uptake  of  certain  low  molecular  weight  dissolved 
organics  only  when  the  microorganism  exists  in  an  advective  flow  field.  As  predicted 
from  a  mass  transfer  model,  fluid  shear  rates  in  natural  systems  are  too  low  to  affect 
bacterial  uptake  of  such  compounds. 
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INTRODUCTION 


The  metabolic  activity  of  heterotrophic  bacteria  attached  to  particles  is  higher  than  that 
of  unattached  bacteria  in  pelagic  ecosystems.  On  a  per-cell  basis,  uptake  by  attached 
bacteria  is  larger  than  uptake  by  free-living  bacteria  for  glucose  (Ibierri  et  al.  1982), 
glucose  and  glutamate  (Kirchman  and  Mitchell  1982),  dissolved  ATP  (Hodson  et  al. 
1981),  phosphate  (Paerl  and  Merkel  1982),  thymidine  (Jeffery  and  Paul  1986),  protein 
hydrosyiate  (Simon  1985),  and  amino  acids  (Bright  and  Fletcher  1983,  Palumbo  et  al. 
1984).  Part  of  this  difference  can  be  attributed  to  the  larger  cell  size  of  attached 
bacteria  compared  with  free-living  bacteria  (Kirchman  1983,  Alldredge  et  al.  1986; 
Simon  1987)  or  to  utilization  of  the  particle  itself.  However,  these  explanations  cannot 
account  for  the  observed  increased  activities  of  bacteria  attached  to  inert  surfaces 
compared  with  similarly-sized  free-living  cells  (Bright  and  Fletcher  1983,  Fletcher  1986). 

The  enhanced  metabolic  activity  of  attached  bacteria  may  also  be  related  to 
differences  in  fluid  environments  experienced  by  attached  versus  free-living  bacteria. 
The  microenvironment  of  free- living  heterotrophic  bacteria  can  be  characterized  as 
laminar  shear.  Unattached  bacteria  must  move  with  the  bulk  fluid,  and  their  cell  size  is 
substantially  smaller  than  the  Kolmogorov  microscale  of  turbulence.  Shear  rates  would 
have  to  approach  10^  s'*  to  reduce  the  microscale  of  turbulence  to  the  size  of  a  1  fim 
bacterium  (Purcell  1978),  which  is  unlikely  in  lakes  and  oceans  where  the  highest  shear 
rates  in  the  upper  few  meters  of  waters  may  only  reach  2-7  s'  *  (Soloviev  et  al.  1988). 
However,  bacteria  attached  to  detritus,  such  as  marine  snow,  can  have  advective  flow 
past  their  surface.  Marine  snow  can  be  highly  porous  and  can  sink  at  velocities 
approaching  200  m  d"  *  (Alldredge  and  Gotschalk  1988).  Theoretical  calculations  (Logan 
and  Hunt  1987,  1988,  Logan  and  Alldredge  1989)  and  experimental  evidence  (Logan 
1987,  Li  and  Ganczarczyk  1987;  Wittier  1986)  indicate  that  velocities  of  fluid  motion 
through  high-porosity  aggregates  approach  the  free-field  velocity.  Therefore,  bacteria 
within  these  aggregates  can  reside  in  an  advective  flow  field  (Logan  and  Hunt  1987). 
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Fluid  motion  should  affect  bacterial  uptake  of  dissolved  organics  analogous  to 
how  fluid  motion  affects  phytoplankton.  Nutrient  uptake  by  phytoplankton  increases 
with  either  fluid  flow  past  the  cell  or  with  fluid  shear  (Caneili  and  Fuhs  1976,  Pasciak 
and  Gavis  1975,  Mierle  1985).  In  laboratory  experiments  with  pure  cultures  of  bacteria, 
Logan  and  Dettmer  (1990)  showed  that  leucine  uptake  by  Zooeloea  ramigera  fixed  in  a 
flow  field  of  1  mm  s_1  was  55-65%  greater  than  uptake  by  suspended  cells.  However, 
shear  rates  as  large  as  50  s“*  did  not  affect  leucine  uptake  (Logan  and  Dettmer  1990). 
The  effect  of  fluid  flow  on  natural  marine  bacterai,  which  differ  from  those  cultivated 
in  the  laboratory,  has  not  been  examined. 

The  purpose  of  this  investigation  was  to  examine  if  fluid  motion  affected  uptake 
by  natural  assemblages  of  marine  bacteria  in  the  same  manner  as  a  pure  culture  and  to 
examine  the  assumptions  of  a  theoretical  model  developed  to  predict  the  effect  of  fluid 
motion  on  uptake.  The  model  is  based  on  the  assumption  that  fluid  motion  enhances 
uptake  by  compressing  the  concentration  boundary  layer  surrounding  the  cell  (Logan  and 
Hunt  1987).  Therefore,  fluid  motion  should  only  increase  uptake  when  uptake  is  not 
already  saturated  and  at  the  maximum  rate.  We  measured  the  uptake  of  two  model 
compounds  (leucine  and  glucose)  by  bacteria  fixed  in  a  uniform  flow  field  and 
suspended  in  a  laminar  shear  field.  We  found  that  bacterial  uptake  of  leucine,  but  not 
glucose,  increased  with  the  rate  of  flow  and  that  the  effect  of  fluid  motion  could  be 
removed  at  elevated  leucine  concentrations.  These  results  imply  that  under  nutrient 
limited  conditions,  attachment  to  detritus  could  increase  bacterial  uptake  if  there  is 
sufficient  advective  flow  past  the  ceil. 

METHODS 

Sampling  sites  and  procedures:  All  samples  were  collected  from  surface  waters  in  the 
Roosevelt  Inlet  on  the  Delaware  Bay  (Lewes,  DE)  during  August,  1989.  Surface  samples 
were  gravity- filtered  through  0.8  Nucleopore  filters  to  remove  bacterivores  and 
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autotrophs.  Samples  were  either  immediately  used  in  experiments  or  kept  overnight 
(aged)  in  the  dark  at  room  temperature. 

Uptake  rates  and  concentrations:  Dissolved  free  amino  acids  (DFAA)  were 
measured  by  high  performance  liquid  chromatography  (HPLC,  Rainin)  using  an 
OPA-reverse  phase  method  (Lindroth  and  Mopper  1979)  with  a  Cl  8-reverse  phase 
column.  The  analytical  procedures  were  essentially  those  of  Fuhrman  and  Bell  (1985) 
with  a-aminobutyric  acid  as  an  internal  standard.  Samples  for  HPLC  analysis  were 
filtered  with  Gelman  Acrodiscs  (25  mm  with  pore  size  of  0.45  pm). 

Leucine  and  glucose  uptake  by  bacteria  was  determined  using  4,5-[^H]leucine  (53 
Ci  mmol-1,  ICN  Laboratories)  and  [^H]glucose  (18  Ci  mmol-1,  ICN  Laboratories).  All 
samples  were  analyzed  using  a  Beckman  3801  Liquid  Scintillation  Counter,  with 
H- number  correction  for  sample  quenching.  Microorganisms  were  enumerated  using  the 
acridine  orange-epifluoresence  direct  count  method  (Hobbie  et  al.  1977). 

The  uptake  rate  of  radiolabeled  glucose  and  leucine  by  suspended  bacteria  in 
undisturbed  fluid  was  determined  prior  to  conducting  experiments  on  the  effects  of 
fluid  motion  on  uptake.  Samples  were  spiked  with  radiolabeled  leucine  or  glucose  at  the 
same  concentrations  used  in  flow  experiments  (-1  or  11  nM).  Samples  (10  ml,  in 
triplicate)  were  withdrawn  at  0,  2  and  10  min,  combined  with  formalin  (2%  final 
concentration)  and  vacuum  filtered  at  250  mm  Hg  through  0.2  pm  nylon  filters  (Poretics 
Corp.).  Filters  were  rinsed  twice  with  5  ml  of  filtered  seawater  and  radioassayed. 

Effect  of  Fluid  Flow:  The  effect  of  fluid  motion  on  bacterial  uptake  was 
examined  by  measuring  uptake  of  leucine  or  glucose  by  cells  that  were  held  on  the 
surface  of  a  filter.  Filter-bound  bacteria  were  prepared  by  filtration  of  the  <0.i~pm 
size  fraction  (10  mi)  through  a  0.2  pm  nylon  filter  (Poretics  Corp.)  supported  by  a  5  pm 
cellulose  acetate  backing  filter  (Millipore  Corp.).  These  filters  were  held  by  stainless 
steel  bases  and  funnels  (25  mm  diameter)  in  a  10-place  vacuum  box  fitted  with  teflon 
valves  (Hoefer  Scientific).  The  nylon  filter  is  a  depth  filter,  and  some  penetration  of 
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cells  into  the  filter  matrix  may  have  occurred.  Based  on  ceil  size  and  filter  area,  we 
calculate  that  ceils  covered  less  than  1%  of  the  filter  surface. 

Flow  past  bacteria  was  obtained  by  pulling  particle-free  seawater  past  the 
filter-bound  bacteria.  Particle-free  media  was  prepared  by  addition  of  radiolabeled 
leucine  or  glucose  to  either  fresh  or  aged  seawater  that  had  been  previously  filtered 
through  a  0.2  Atm  polycarbonate  filter  (Poretics  Corp.).  This  solution  was  applied  to  the 
filter-bound  bacteria  with  a  syringe  with  a  0.2  nm  in-line  filter  (Gelman  Corp.). 
Different  fluid  velocities  were  obtained  by  varying  the  vacuum;  incubation  times  were 
kept  between  1-2  minutes  by  varying  flow  volumes.  After  all  the  fluid  had  been  pulled 
through  the  filter,  bacteria  were  immediately  killed  using  2  ml  of  2%  formalin,  and 
rinsed  with  an  additional  10  mi  of  0.2-Aim  filtered  seawater.  Filters  were  added  to  7  ml 
of  scintillation  cocktail  (Aquasol  II)  and  radioassayed. 

Effect  of  Fluid  Shear  A  defined  laminar  shear  environment  was  obtained  using  a 
rotating  cylinder  apparatus  (van  Duuren  1968).  Our  device  consisted  of  two  vertically 
aligned  concentric  acrylic  cylinders  separated  by  a  9  mm  gap.  The  fixed  inner  cylinder 
was  8  cm  in  diameter  and  15  cm  high.  The  outer  cylinder  rotated  at  a  fixed  speed 
producing  a  uniform  shear  rate  between  the  cylinders. 

Radiolabeled  substrates  were  added  to  500  ml  of  the  aged  bacterial  size  fraction 
(<0.8  Aim),  and  400  ml  of  this  solution  was  poured  into  the  rotating  cylinder  annulus. 
The  cylinder  was  started,  and  time  zero  samples  were  obtained  by  withdrawing  three 
10-ml  samples  from  the  remaining  undisturbed  100  ml.  After  10  min  the  rotating 
cylinder  was  stopped  and  the  inner  cylinder  withdrawn.  Three  10-ml  samples  were 
withdrawn  from  the  rotating  cylinder  and  from  the  undisturbed  bacterial  size  fraction. 
Formalin  (2%  final  concentration)  was  added  and  the  sample  was  filtered  through  0.2  fim 
nylon  filters.  Uptake  was  calculated  as  the  difference  between  initial  and  final 
radioactivity  of  the  filters. 

Arsenate  inhibition  of  leucine  and  glucose  uptake:  To  begin  to  determine 
whether  or  not  leucine  and  glucose  are  transported  by  similar  mechanisms,  the  effect  of 
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arsenate  on  uptake  was  examined.  The  monopotassium  salt  of  arsenic  acid  (Sigma)  was 
added  directly  to  an  aged  <0.8  ^m  size  fraction  to  avoid  leucine  or  glucose 
contamination.  After  addition,  the  seawater  was  mixed  and  incubated  for  3  min. 
[^H]leucine  or  [^H]glucose  was  then  added  and  uptake  followed  (in  duplicate)  over  time. 
Filters  were  rinsed  once  with  4  ml  of  filtered  seawater  and  radioassayed.  Controls  did 
not  receive  any  arsenate. 

Kinetic  model:  The  effect  of  fluid  motion  on  microbial  kinetics  can  be 
determined  using  a  mass  transfer  analysis  which  is  summarized  below.  Details  of  the 
model  are  described  elsewhere  (Logan  and  Alldredge  1989,  Logan  and  Dettmer  1990). 
Assuming  substrate  uptake  is  proportional  to  substrate  concentration,  C,  the  rate  of 
nutrient  utilization  can  be  calculated  as: 

-1  dC 

- =,  k4Tra2(C-C.)  (1) 

N  dt 

where  N  is  the  cell  concentration,  k  is  a  mass  transfer  coefficient,  a  is  the  cell  radius, 
and  C  and  Cs  are  the  bulk  and  cell-surface  substrate  concentrations,  respectively. 
According  to  the  Michaelis-Menten  description  of  enzyme  kinetics,  a  molecule  reversibly 
and  repeatedly  binds  to  an  enzyme  before  successful  product  production.  Since  few 
enzymes  achieve  "kinetic  perfection"  (Stryer  1981),  or  operate  at  the  limit  of 
diffusion-controlled  uptake  rates,  the  substrate  concentration  at  the  cell  surface  is: 

Cs  *  (1  -  a)  C  (2) 

where  a,  the  collision  efficiency,  is  limited  to  the  range  of  0£a£l.  The  mass  transfer 
coefficient  can  be  reduced  in  proportion  to  enzyme  coverage  as  long  as  the  fraction  of 
the  surface  active  in  transport,  f,  is  small  (Dudukovic  and  Mills  1985).  By  defining 
Eg-af,  where  0£Eg£l,  we  obtain: 

-1  dC 

- -  4tt  a  D  Sh  EB  C  (3) 

N  dt 

where  Sh  is  the  dimensionless  Sherwood  number,  defined  as  Sh-ka/D,  and  D  is  the 
substrate  diffusivity. 
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For  a  spherical  cell  in  stagnant  fluid,  the  Sherwood  number  is  unity.  The 
Sherwood  number  has  no  physical  limits,  but  high  fluid  velocities  past  the  cell  will 
ultimately  damage  the  cell  and  decrease  nutrient  uptake.  The  Sherwood  number 
correlation  used  in  this  study  to  predict  the  effect  of  fluid  flow  past  cells  (Logan  and 
Alldredge  1989)  was: 

Shu  »  1  +  45  Re0-57  (4) 

where  Re«au/v  is  the  cell  Reynolds  number,  u  is  the  fluid  velocity  past  the  cell,  and  v 
is  the  kinematic  fluid  viscosity.  For  a  cell  in  laminar  shear  flow,  Frankel  and  Acrivos 
(1968)  proposed: 

ShG  -  1  +  0.26  Pe0-5  (5) 

where  Pe  is  the  Peclet  number,  defined  as  Pe-a^G/D-ccl,  and  G  is  the  shear  rate. 

RESULTS 

Effects  of  fluid  flow  on  [^HJleucine  uptake:  Fluid  motion  is  predicted  to  increase 
microbial  uptake  only  if  transport  is  not  saturated.  During  our  experiments,  leucine 
transport  was  probably  saturated  in  surface  waters  of  Roosevelt  Inlet  since  the  leucine 
concentration  was  45  nM  (August  7)  and  Ks  for  leucine  uptake  is  probably  <1  nM.  To 
reduce  leucine  concentrations  and  obtain  non-saturated  uptake  kinetics,  we  aged  the 
bacterial  size  fraction  (<0.8  fim  filtrate)  for  24  hr  in  the  dark;  in  the  absence  of  leucine 
input  from  other  microorganisms,  bacterial  uptake  reduced  leucine  concentrations  to  i.7 
nM.  Leucine  concentrations  in  particle-free  seawater  were  near  the  detection  limit  (<1 
nM),  indicating  that  the  filtered  water  used  for  the  flow  experiments  was  not 
contaminated. 

At  low  leucine  concentrations,  fluid  flow  past  cells  increased  [^Hjleucine  uptake 
by  a  factor  of  8  compared  to  uptake  at  low  fluid  velocities  (Fig.  1).  [^H]Leucine 
uptake  increased  from  0.65  pg  ml-1  min"1  at  3  ra  d"1  to  a  maximum  of  5.4  pg  ml"1 
min"1  at  -20  m  d"1.  At  higher  fluid  velocities,  [^H]Ieucine  uptake  decreased  to  an 
average  of  4.5  pg  ml"1  min"1  for  fluid  velocities  between  20-70  m  d"1.  The  maximum 
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rate  of  leucine  uptake  in  flow  experiments  was  63%  less  than  the  average  uptake  of 
suspended  cells  (14.5±0.7  pg  ml"1  min"1,  ±SD).  Reduced  uptake  by  cells  on  the  filter 
was  observed  in  other  experiments.  The  reasons  for  these  decreases  are  discussed  below. 

To  test  our  assumption  that  enhanced  activity  required  unsaturated  uptake 
kinetics,  we  repeated  the  flow  experiment,  but  added  radiolabeled  leucine  at  a 
concentration  expected  to  saturate  uptake  kinetics.  The  concentration  of  leucine  in  aged 
samples  for  this  experiment  was  1.1  nM,  and  was  below  detection  (<1  nM)  in 
bacteria-free  samples.  When  10.8  nM  [^H]leucine  was  added  to  the  bacteria-free  sample 
used  in  the  flow  experiments,  no  increase  in  uptake  was  observed  with  fluid  motion 
(Fig.  2).  [3H]Leucine  uptake  by  suspended  bacteria  averaged  18.5±2.5  pg  ml"1  min"1. 
Filter  bound  bacteria  in  flow  fields  ranging  from  8-50  m  d"1  had  an  average  of  13.5±1.0 
pg  ml"1  min"1,  or  27%  less  than  suspended  samples.  This  indicates  that  flow  increases 
up  take  only  when  uptake  is  not  saturated. 

Effects  of  fluid  flow  on  [^Hlglucose  uptake:  We  also  examined  the  effect  of 
fluid  motion  on  [3H]gIucose  uptake  using  aged  samples  from  the  Roosevelt  Inlet.  We 
were  unable  to  detect  any  increase  in  fluid  motion  on  uptake  (Fig.  3).  [3H]glucose 
uptake  by  suspended  averaged  3.9±0.3  pg  ml"1  min'1,  when  glucose  was  added  at  a 
concentration  of  1  nM.  At  fluid  velocities  between  8-48  m  d"1,  glucose  uptake  was  an 
average  of  3.1±0.4  pg  ml"1  min"1,  or  about  23%  lower  than  suspended  samples. 

Effect  of  fluid  shear  on  (^Hlleucine  uptake:  Fig.  4  shows  the  ratio  of  leucine 
uptake  by  suspended  ceils  to  uptake  by  cells  suspended  in  a  homogeneous  laminar  shear 
field.  At  shear  rates  of  0-2.1  s'1,  fluid  shear  did  not  appreciably  affect  leucine  uptake 
of  either  fresh  or  aged  samples  by  more  than  ±10%.  These  results  are  consistent  with 
mass  transfer  calculations  using  Eq.  3  and  5.  At  a  shear  rate  of  2  s'1,  the  Sherwood 
number  is  1.014  for  a  ceil  0.4  n m  in  radius.  The  data  shown  in  Fig.  4  are  not 
sufficiently  precise  to  detect  a  1.4%  increase  ic  leucine  uptake,  and  indicate  that  leucine 
uptake  may  actually  have  decreased  by  up  to  10%  in  sheared  samples. 
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Arsenate  inhibition:  To  determine  whether  the  uptake  systems  for  leucine  and 
glucose  were  similar,  we  examined  the  effect  of  arsenate  on  leucine  and  glucose  uptake, 
by  suspended  bacteria.  Arsenate  inhibits  ATPase  which  hydrolyzes  ATP  to  provide 
energy  for  some  transport  systems  (Ames  1986).  Arsenate  (10  mM)  completely  inhibited 
[^H]leucine  transport  over  5  min  incubations  (Fig  5a).  In  contrast,  uptake  of  [^H]glucose 
was  inhibited  by  only  50%  (Fig  5b).  These  results  suggest  that  the  mechanism  for 
leucine  and  glucose  uptake  differ. 

DISCUSSION 

Fluid  motion  can  substantially  increase  uptake  by  microorganisms,  but  this  effect  is  a 
function  of  several  factors,  including  the  substrate  type  and  its  concentration  relative  to 
uptake  kinetics.  Marine  bacteria  held  on  filters  in  a  flow  stream  of  10  m  d-1  assimilated 
leucine  8  times  faster  than  bacteria  under  similar  conditions  in  the  absence  of  fluid 
motion.  This  enhancement  was  only  observed  when  leucine  uptake  was  not  saturated. 
In  order  for  fluid  motion  to  increase  uptake,  concentration  gradients  must  be  present  at 
the  cell  surface.  At  high  leucine  concentrations  (-11  nM),  fluid  motion  had  no  effect  on 
leucine  uptake  probably  because  of  the  lack  of  a  concentration  boundary  layer. 

The  observed  8-fold  increase  in  [^H]leucine  uptake  at  velocities  between  1-20  m 
d-1  is  much  larger  than  expected  from  the  mass  transfer  analysis  using  Eq.  4.  At  a 
velocity  of  20  m  d*1,  the  Sherwood  number  for  a  cell  0.4  fim  in  radius  is  only  20% 
larger  than  for  a  suspended  cell,  indicating  that  fluid  flow  should  only  increase  uptake 
by  a  factor  of  1.2.  The  8-fold  increase  is  not  only  larger  than  predicted  by  a  mass 
transfer  analysis,  it  is  also  inconsistent  with  other  experimental  results.  Pure  cultures  of 
bacteria  (Zoogloea  ramigera)  fixed  in  a  flow  field  of  1  mm  s-1  increased  by  a  factor  of 
1.55-1.65  compared  to  suspended  cultures  (Logan  and  Dettmer  1990).  Canelli  and  Fuhs 
(1976)  observed  a  2.6-fold  increase  in  phosphorus  uptake  by  the  diatom  Thalassiosira 
fluviatilis  (7  /im  radius)  held  on  a  filter  surface  in  a  velocity  field  of  35  m  d“*.  Since 
the  Sherwood  number  is  proportional  to  the  ceil  radius,  we  expected  the  magnitude  of 
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the  increase  in  uptake  by  bacteria  would  be  smaller,  not  larger,  than  increases  observed 
for  diatoms  by  Canelli  and  Fuhs. 

Bacterial  uptake  rates  of  leucine  and  glucose  by  suspended  bacteria  were 
consistently  larger  than  uptake  rates  determined  for  filter-bound  bacteria.  In  other 
experiments  using  pure  cultures  of  microbes  (Logan  and  Dettmer  1990),  there  was  no 
difference  between  uptake  rates  of  suspended  cells  and  filter-bound  cells  held  in  fluid 
flow  fields  at  low  velocities.  It  is  likely  that  natural  bacteria,  which  are  smaller  than 
pure  cultured  bacteria  previously  examined  in  laboratory  experiments,  became  partially 
enmeshed  in  the  filter  material,  reducing  the  exposed  surface  area  of  the  cell.  In  the 
flow  experiments,  uptake  rates  at  higher  fluid  velocities  (20-70  m  d'*)  averaged  ~65% 
less  than  suspended  experiments.  This  indicates  that  as  much  as  65%  of  the  surface  area 
of  the  cell  used  for  leucine  uptake  was  covered  by  the  filter  during  flow  experiments. 
Reduced  cell  surface  area  by  the  filter,  therefore,  could  explain  the  observed  reduction 
in  leucine  uptake  by  cells  held  on  filters. 

The  area  covered  by  outer  membrane  proteins  involved  in  leucine  transport 
(leucine  porters)  is  important  for  examining  the  impact  of  flow  on  leucine  utilization  and 
for  understanding  how  these  bacteria  survive  in  the  low  o.ganic  v.wucwuudtions  of 
seawater.  This  surface  area  can  be  calculated  using  Eq.  3.  For  5.3  x  106  cell  ml-1,  and 
assuming  a  leucine  diffusivity  of  7  x  10"**  cm^  s'*  (Logan  and  Dettmer  1990)  and  a  cell 
radius  of  0.4  ^im.  Eg  equals  0.09  for  an  uptake  rate  of  14.5  pg  ml"*  min'*  at  a  leucine 
concentration  of  1.1  nM.  If  we  assume  the  porter  is  100%  efficient  (a-1),  this  Eg 
indicates  that  9%  of  the  surface  area  of  free- living  cells  is  covered  by  leucine  porters. 
This  percent  coverage  is  close  to  the  surface  area  of  porins,  which  are  abundant 
nonspecific  pores  or  channels  found  in  Gram- negative  bacteria  (e.g.  £.  coli).  Assuming 
that  each  £.  coli  cell  contains  approximately  10^  porins  with  a  surface  area  of  0.9  nm?- 
(Nikaido  and  Vaara  1987)  and  a  total  cell  surface  area  of  6.7  fi m^,  approximately  13%  of 
the  cell  surface  is  occupied  by  porins.  The  surface  area  of  specific  leucine  porters  in  £. 
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coli  is  probably  much  less  than  for  marine  bacteria,  although  we  are  unaware  of  any 
published  estimates. 

We  did  not  observe  any  effect  of  fluid  motion  on  uptake  of  glucose  uptake  in 
either  new  or  aged  water  samples.  Since  glucose  concentrations  were  not  measured,  we 
cannot  rule  out  the  possibility  that  glucose  uptake  was  saturated  and  that  flow  actually 
affects  glucose  uptake  like  leucine.  However,  saturation  of  glucose  uptake  seems 
unlikely,  especially  in  aged  samples  since  concentrations  of  labile  compounds  like  amino 
acids  are  substantially  reduced  in  24  h  in  the  bacterial  size  fraction.  Alternatively,  we 
hypothesize  that  the  different  effects  of  flow  on  uptake  of  glucose  and  leucine  is  due  to 
different  transport  mechanisms.  The  arsenate  experiment  provides  some  support  for  this 
argument.  Amino  acid  transport  is  often  mediated  by  periplasmic  binding  proteins 
(Furlong  1987)  fueled  by  ATP  hydrolysis  (Ames  1986),  and  thus  dependent  on  ATPase. 
Arsenate,  a  ATPase  inhibitor,  should  block  transport  of  compounds  such  as  leucine 
(Ames  1986).  That  arsenate  did  not  completely  inhibit  glucose  uptake  suggests  that  a 
binding  protein  is  not  involved  in  glucose  uptake.  Other  mechanisms  for  glucose 
transport  do  not  involve  a  specific  outer  membrane  protein  (e.g.,  the  PTS  system; 
Meadows  et  al.  1987,  Delong  and  Yayanos  1987).  According  to  these  mechanisms, 
glucose  enters  the  periplasmic  space  of  Gram-negative  bacteria  (the  most  common 
bacteria  in  marine  systems)  via  a  nonspecific  porin.  Passive  diffusion  through  porins 
requires  a  concentration  gradient  across  the  membrane.  This  gradient,  and  not  a 
concentration  boundary  layer  outside  the  cell,  could  provide  the  largest  resistance  to 
mass  transfer  into  the  cell.  Thus,  flow  could  enhance  uptake  mediated  by  specific 
binding  proteins  by  reducing  the  concentration  boundary  layers  surrounding  a  ceil,  but 
not  uptake  where  flux  is  limited  by  a  gradient  across  the  outer  cell  membrane. 

Our  uptake  experiments  conducted  on  marine  bacteria  at  shear  rates  below  2  s'1 
indicate  that  fluid  shear  is  ineffective  at  increasing  microbial  uptake  of  small  molecular 
weight  compounds  such  as  leucine.  Since  free-living  marine  bacteria  are  much  smaller 
than  the  Kolmogorov  microscale  of  turbulence,  fluid  shear  is  ineffective  in  reducing 
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concentration  gradients  at  the  ceil  surface.  Our  observations  are  supported  by  mass 
transfer  correlations  (Eq.  5),  which  predict  a  maximum  increase  of  1.4%  at  shear  rates 
below  2  s"1  and  are  consistent  with  previous  theoretical  predictions  (Purcell  1978,  Logan 
and  Hunt  1987). 

In  conclusion,  wc  demonstrated  that  fluid  flow  of  only  20  m  d-1  may  be 
sufficient  to  cause  a  several-fold  increase  in  uptake  rates  of  selected  compounds  by 
attached  bacteria.  Since  sinking  velocities  of  particles  may  reach  200  m  d_1  (Alldredge 
and  Gotschalk  1988),  attached  bacteria  may  have  a  selective  advantage  over  free-living 
bacteria  in  oligotrophic  environments.  Increased  mass  transfer  to  attached  cells  may 
partially  explain  why  uptake  by  attached  bacteria  in  natural  environments  is  higher  than 
that  of  free-living  bacteria. 
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Abstract 

Blooms  of  chain-forming  diatoms  commonly  flocculate  into 
centimeter-sized  aggregates  of  living,  vegetative  cells  follow¬ 
ing  nutrient  depletion  in  surface  waters  off  southern  Califor¬ 
nia.  We  examined  the  hypothesis  that  diatom  cells  within 
aggregates  experience  increased  nutrient  uptake  relative  to 
unattached  cells.  We  measured  in  situ  settling  velocities  of 
49  to  190  md’1  and  calculated  porosities  of  0.99931  to 
0.99984  ( ±  <  0.03%)  for  12,  newly-formed  diatom  floes 
ranging  from  0.19  to  4.2  cm3  in  volume  and  7  to  22  mm  in 
equivalent  spherical  diameter.  Using  permeability -porosity 
relationships,  we  calculated  intra-aggregate  flow  velocities 
of  20  to  160 yum  s'1.  Although  subject  to  considerable  un¬ 
certainty,  a  Relative  Uptake  Factor  analysis  based  on  mass 
transfer  equations  suggests  that  diatoms  fixed  within 
aggregates  undergoing  gravitational  settling  can  take  up  nu¬ 
trients  up  to  2.1  ±  0.4  times  faster  than  unattached  diatoms 
experiencing  laminar  shear.  Increased  nutrient  uptake  by 
aggregated  diatoms  may  be  important  in  understanding  the 
reasons  for  diatom  floe  formation. 


Introduction 

Blooms  of  chain-forming  diatoms  often  disappear  rapidly 
from  surface  ocean  waters  following  nutrient  depletion 
(Guillard  and  Kilham  1977,  Smetacek  1985).  This  mass  sedi¬ 
mentation  occurs  via  flocculation  of  chains  into  loose, 
centimeter-sized  aggregates  of  living,  vegetative  cells  (All¬ 
dredge  and  Silver  1988,  Kranck  and  Milligan  1988)  which 
sink  at  rates  of  50  to  200  md'1,  one  to  two  orders  of  magni¬ 
tude  higher  than  unaggregated  chains  (Alldredge  and  Got- 
schalk  in  press).  Mass  sinking  of  diatoms  has  been  equated 
with  mass  mortality  (Walsh  1983),  and,  at  first  glance,  mech¬ 
anisms  resulting  in  the  rapid  settlement  of  viable  cells  out  of 
the  euphotic  zone  appear  maladaptive.  However,  Smetacek 
(1985)  suggested  that  mass  flocculation  removes  diatoms 
from  warm,  nutrient-depleted  surface  waters,  where  their 
viability  is  low  into  dark,  colder  waters  where  they  may 


remain  viable  for  long  periods  of  time  as  resting  spores  and 
vegetative  resting  stages.  Later  resuspension  of  these  resting 
stages  back  into  the  euphotic  zone  would  reseed  surface 
waters  with  diatoms  when  conditions  were  again  favorable 
for  population  growth. 

While  this  evolutionary  explanation  of  diatom  floccula¬ 
tion  is  supported  by  existing  literature  on  diatom  biology,  it 
does  not  preclude  a  more  immediate  adaptive  advantage  for 
mass  flocculation.  There  are  several  possible  immediate  ad¬ 
vantages  for  cell  growth  within  the  microhabitat  of  a  porous 
aggregate,  including  interactions  between  adjacent  micro¬ 
organisms.  (i.e..  commensalism,  mutualism,  and  exchange 
of  genetic  material)  and  potential  protection  from  certain 
predators.  The  most  important  advantage  in  highly  potous 
aggregates,  however,  may  be  increased  nutrient  uptake  (Lo¬ 
gan  and  Hunt  1987).  Many  types  of  cells,  including  yeasts 
and  bacteria,  aggregate  when  nutrients  stressed  (see  Calleja 
1984)  and  diatom  flocculation  in  nature  also  appears  to  be 
triggered  by  nutrient  limitation  (Smetacek  1985).  Moreover, 
cultures  of  Chaetoceros  debilis,  a  representative  of  the  most 
common  genus  in  natural  diatom  fiocs,  flocculate  under 
gentle  rotation  in  the  laboratory  only  when  they  reach  sta¬ 
tionary  growth  phase,  (i.e.,  nutrient  limitation)  and  not 
before  (Alldredge  unpublished  data).  Logan  and  Hunt 
(1987)  predicted  that  advective  flow  through  porous  bacte¬ 
rial  aggregates  would  increase  the  overall  uptake  of  large 
molecular  weight,  dissolved  organic  molecules  by  attached, 
relative  to  unattached,  bacteria.  Moreover,  Munk  and  Riley 
(1952)  suggested  that  phytoplankton  cells  experiencing 
faster  sinking  velocities  also  have  accelerated  nutrient  ab¬ 
sorption  and  increased  growth  rates  despite  decreased  con¬ 
centrations  of  nutrients.  Although  Munk  and  Riley  (1952) 
did  not  consider  the  case  of  phytoplankton  in  aggregates, 
recent  experiments  on  the  absorption  of  a  large  molecular 
weight  dye  by  different  types  of  marine  snow  support  the 
existence  of  bulk  fluid  motion  through  natural  marine 
aggregates  (Logan  1987)  and  suggests  that  diatoms  in  sink¬ 
ing  aggregates  experience  relatively  more  rapid  riuid  flow 
than  unattached  diatoms. 
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In  this  study  we  examined  the  hypothesis  that  diatom 
cells  within  aggregates  experience  increased  overall  trans¬ 
port  ol'  nutrients  relative  to  unattached  cells,  i.e..  that  an 
immediate  adaptive  advantage  of  diatom  flocculation  is  in¬ 
creased  nutrient  uptake  by  the  aggregated  cells.  We  used 
measured  sinking  velocities  and  porosities  to  calculate  intra- 
aggiCuait  flow  velocities  p«oi  cells  fixed  wifi.m  natural  di¬ 
atom  floes.  We  then  used  mass  transfer  theory  to  predict  the 
overall  nutrient  uptake  of  cells  within  natural  diatom  floes 
relative  to  unattached  cells. 


Methods 

Experimental 

The  settling  velocities  of  12  undisturbed,  newly-formed  di¬ 
atom  floes  were  measured  directly  in  situ  by  SCUBA  divers 
in  surface  waters  of  the  Santa  Barbara  Channel.  California 
(34  23'N.  119  50'W)  in  July  1986.  Settling  velocities  were 
determined,  as  described  in  Alldredge  and  Gotschalk  (1988). 
by  measuring  the  time  required  for  each  aggregate  to  sink  to 
a  spot  of  neutrally  buoyant  fluorescein  dye  placed  3  cm 
below  the  aggregate  in  the  water  column.  Aggregate  volume 
was  determined  from  photographs  taken  in  suu  and  each 
aggregate  was  collected  for  gravimetric  determination  of  dry 
weight. 

Porosity,  p.  was  calculated  directly  from  measurements 
of  aggregate  volume  and  dry  weight  using  the  equation 


thng.  Masliyah  and  Polikar  ( 1980)  determined  the  empirical 
expression  of  the  drag  coefficient  as: 

24  Q 

C0  =  -  -  [1  +  0.0853  Re"  093  0  1,15  “’].  (3) 

Re 

where  Re  —  L'pd  t.  i  is  the  kinematic  viscosity  of  the  fluid, 
u  =  log10(Re).  and  Q  is  determined  from 

2c2  [1  -  tan  hi  cl 

Q  =  — -c - — (4) 

2c2  +  3[1  -  tan  he  c] 

where  c  is  the  dimensionless  size  of  the  aggregate,  defined  as 
l  =  o,  K  1  2.  af  is  the  aggregate  radius,  and  K  is  the  permea¬ 
bility  of  the  aggregate 

The  permeability  of  highlv-porous  aggregates  was  de¬ 
fined  as  a  function  of  the  aggregate  structure  and  porosity 
using  (Davies  1952): 

K~ 1  =  16u(._"  °ll  +  5603).  (5) 

where  ac  is  the  cell  radius,  and  <i>  is  the  occupied  volume 
fraction  of  the  aggregate,  equal  to  the  quantity  ( 1  —  p\.  We 
used  the  Davies  correlation  since  this  expression  was  used  by 
Masilyah  and  Polikar  (1980)  in  settling  experiments  to  vali¬ 
date  the  existence  of  flow  through  permeable  aggregates 
with  0.2  <  Re  <  120. 

The  average  intra-aggregate  volocity,  u,  through  a  per¬ 
meable  spherical  aggregate  undergoing  gravitational  settling 
(Adler  1981)  is: 


P=  1  - 


u  =  V 

a, 


where  W  is  the  dry  weight  of  the  aggregate.  V  is  aggregate 
volume  and  P  is  the  density  of  the  diatoms  composing  the 
aggregate,  here  assumed  to  be  1.1047  gem  3.  This  value  was 
calculated  from  ambient  seawater  density  (1.0247  gem'3. 
Alldredge  and  Gotschalk  1988)  assuming  an  excess  density 
for  nutrient-depleted  living  diatom  cells  of  0.08  gem'3 
( Eppley  et  al.  1 967).  Species  composition  and  mean  cell  sizes 
of  diatoms  in  floes  were  determined  by  counting  and  mea¬ 
suring  cells  in  eight  additional  floes  collected  at  the  same 
time. 


Theoretical 

Measured  sinking  velocities  were  compared  to  settling  veloc¬ 
ities  of  highly  porous  spherical  aggregates,  determined  from 
a  force  balance  as 


l 


■2 

P 


8  q(P-Pf)d 
6CDPr 


(2) 


where  Up  is  the  predicted  settling  velocity,  q  is  the  accelera¬ 
tion  due  to  gravity,  P  is  the  particle  density,  Pf  is  the  fluid 
density,  d  is  the  aggregate  diameter,  and  CD  is  the  drag 
coefficient.  For  permeable  aggregates  with  Reynolds  num¬ 
bers  (Re)  between  7  and  120  undergoing  gravitational  set- 


where  U  is  the  measured  in  situ  sinking  velocity,  and  a*  is  a 
function  of  the  dimensionless  variable  c  contained  in  Adler 
(1981).  When  a  highly  porous  aggregate  is  permeable, 
streamlines  cross  the  surface  of  the  aggregate,  and  advect 
fluid  past  diatoms  fixed  within  the  interior  of  the  aggregate. 
Therefore,  a  diatom  within  a  rapidly  sinking  permeable 
aggregate  potentially  experiences  higher  fluid  flow  past  its 
surface  than  does  an  unattached  cell,  which  either  sinks  very 
slowly  or  moves  entirely  with  the  bulk  fluid. 

The  rate  of  nutrient  absorption  by  a  cell  is  a  function  of 
nutrient  concentration  and  flow  fieldifluid  flow  past  a  cell 
increases  the  rate  of  nutrient  absorption  by  the  cell  (Munk 
and  Riley  1952).  Following  the  analysis  of  Logan  and  Hunt 
(1987),  we  incorporated  this  hydrodynamic  effect  into  a  first- 
order  kinetic  model  using  a  mass  transfer  coefficient  to  in¬ 
crease  cell  kinetics  with  increased  fluid  flow  past  the  surface 
of  the  cell.  Since  attached  cells  within  aggregates  are  exposed 
to  lower  nutrient  concentrations  than  unattached  cells,  nu¬ 
trient  removal  by  all  cells  within  the  aggregate  was  com¬ 
pared  to  nutrient  removal  by  an  equivalent  number  of  cells 
freely  dispersed  in  the  bulk  fluid  using  a  Relative  Uptake 
Factor,  as  described  below. 

The  rate  of  mass  transfer  to  a  single  cell,  Pc ,  in  any  fluid 
environment  (Logan  and  Hunt  1987)  is: 


Pc  =  4nacD  EBShC 


(7) 
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where  D  is  the  nutrient  diffusivity.  assumed  to  be 
10  'em  2  s  '.  and  C  is  the  nutrient  concentration.  The 
Sherwood  number  is  a  dimensionless  ratio  of  the  rate  of 
mass  transfer  by  advection  to  the  rate  of  mass  transfer  by 
diffusion,  and  is  defined  here  as  Sh  =  k u  D.  where  k  is  the 
mass  transfer  coefficient.  For  a  spherical  object  of  radius  uc. 
‘he  Slier .vood  number  is  unit)  for  mass  transfer  bv  diffusion 
alone,  and  greater  than  unity  for  increased  mass  transfer  due 
to  fluid  motion  around  the  cell.  The  kinetics  of  nutrient 
uptake  are  incorporated  in  the  transport  model  using  a  col¬ 
lector  efficiency.  EB.  which  defines  uptake  as  that  fraction  of 
collisions  between  substrate  and  cells  that  are  successful;  for 
instantaneous  reactions.  E„  =  1.  and  in  the  absence  of  nu¬ 
trient  uptake.  Eg  =  0.  The  mass  transfer  calculations  pre¬ 
sented  below  use  EB  =  0.025.  as  determined  by  Logan  ( 1986) 
from  experimental  data  reported  by  Canelli  and  Fuhs  (1976). 

In  order  to  examine  mass  transfer  to  attached  and  un¬ 
attached  diatoms,  two  Sherwood  numbers  are  required:  one 
for  diatoms  fixed  in  a  flow  field,  and  a  second  for  unattached 
diatoms  experiencing  laminar  shear.  Sherwood  numbers  for 
diatoms  within  a  permeable  aggregate,  or  for  unattached 
diatoms  sinking  in  undisturbed  fluid,  were  derived  from  data 
reported  by  Canelli  and  Fuhs  (1976)  on  phosphorus  uptake 
by  cells  of  Thalassiosiru  Jluviatilis  fixed  in  a  flow  field.  Using 
their  data,  we  empirically  determined  the  following  correla¬ 
tion  for  the  Sherwood  number  as  a  function  of  the  Reynolds 
number: 

Sh  =  1  +  45  Re0  5" .  (8) 

Many  mass  transfer  relationships  indicate  the  Shervood 
number  is  proportional  to  the  Reynolds  number  to  the  0.5 
power,  and  our  correlation  is  in  good  agreement  with  this 
value,  since  the  standard  error  of  the  power  in  Eq.  8  s  0.08. 
In  Fig.  1  we  compare  the  phosphorus  flux  data  reported  by 
Canelli  and  Fuhs  (their  experiments  2  and  3.  in  Fig.  2)  with 
fluxes  predicted  using  the  above  correlation  (Eq.  8).  At  lower 
concentrations  of  phosphorus  ( 1 0  /tg  I  *  1 ),  flux  values  are  less 
than  10” 11  /tg-P//m  2  min  ',  and  the  observed  and  predict¬ 
ed  fluxes  are  in  relatively  close  agreement.  However,  at 
larger  phosphorus  concentrations  (100/rgl  l),  fluxes  are 


Fig.  1.  Phosphorus  fluxes  (IO~l2//g-P  //m'2min'')  calculated 
using  Eq.  8  (continous  line)  versus  data  reported  by  Canelli  and 
Fuhs  (1976):  ■  -  their  experiment  2  (ac  =  7.3  /rm.  C  =  10  //g-P  1"  '); 
•  their  experiment  3  (a,  =  8.6  ^m.  C  =  100  /(g-P  1' ') 


underestimated  for  fluid  velocities  less  than  200 /im  s  '.  and 
overestimated  for  velocities  greater  than  200 /tm  s  1 .  Canelli 
and  Fuhs'  data  suggest  that  a  saturation  value  was  reached 
in  their  experiments. 

The  fluid  environment  of  an  unattached  diatom  can  be 
characterized  as  laminar  shear,  since  diatoms  freely  sus¬ 
pended  in  the  ocean  are  much  smaller  than  the  Kolmogorov 
microscale  of  turbulence.  The  Sherwood  number  for  a 
sphere  in  laminar  shear  flow.  Sh(i.  (Frankel  and  Acrivos 
1968)  is: 

Sh(;  =  1  +  0.26  Pe^  2  (9) 

where  Pe  is  the  dimensionless  Peelet  number,  defined  as 
Pe,;  =  a2  G  D .  and  G  is  the  fluid  shear  rate.  Eq.  9  is  valid  for 
Peelet  numbers  much  less  than  unity,  and  is  only  used  for 
calculations  of  the  rate  of  mass  transfer  to  unattached  di¬ 
atoms. 

The  ratio  of  the  overall  nutrient  uptake  by  aggregated 
cells  to  the  nutrient  uptake  by  an  equivalent  number  of 
unattached  cells  is  defined  as  the  Relative  Uptake  Factor 
(Logan  and  Hunt  1987).  For  nutrient  uptake  by  cells  in  a 
permeable  aggregate  undergoing  gravitational  settling,  com¬ 
pared  to  unattached  cells  suspended  in  laminar  shear  flow, 
the  Relative  Uptake  Factor.  va,  is: 

Sh„  [1  -exp( -</>„)] 

/a  Sh0  4>a 

where  <j>a.  the  dimensionless  advective  Thiele  Modulus,  is 


2 ‘'A, 

u 


(ID 


and  af  is  the  aggregate  radius,  and  ku  is  the  first  order  uptake 
rate  constant  describing  nutrient  removal  within  the 
aggregate.  From  the  rate  of  uptake  per  cell  (Eq.  7),  ku  is 

3E,Sh.D(l-£),  (12) 

a- 

The  ratio  of  Sherwood  numbers  in  Eq.  10  accounts  for  the 
different  fluid  mechanical  environment  of  the  unattached 
cells  (laminar  shear)  and  the  cells  within  the  aggregate  (fixed 
in  a  flow  field).  The  Thiele  modulus,  as  defined  in  Eq.  1 1,  is 
a  ratio  of  the  rate  of  mass  transfer  to  attached  diatoms,  to 
the  rate  of  fluid  flow.  Since  nutrients  are  removed  by  cells 
as  fluid  flows  through  the  aggregate,  the  Thiele  modulus 
accounts  for  decreased  nutrient  concentrations  in  the 
aggregate  interior. 

For  growth  within  the  aggregate  to  be  advantageous  to 
all  cells  within  the  aggregate,  the  Relative  Uptake  Factor 
must  be  greater  than  unity.  For  v0  <  1,  only  cells  near  the 
surface  of  the  aggregate  would  benefit  from  advective  flow 
through  the  aggregate.  Since  first  order  kinetics  are  implied 
for  all  substrate  concentrations  below  the  bulk  nutrient  con¬ 
centration.  a  Relative  Uptake  Factor  analysis  does  not  re¬ 
quire  that  the  bulk  nutrient  concentration  be  known  since 
uptake  is  measured  relative  to  unattached  cells  in  the  bulk 
fluid. 
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Results 

Taxonomic  composition  offices  expressed  as  percent  of  total 
cell  number  is  shown  in  Table  1.  Three  species  of  diatoms. 
Chaetocerns  radiums.  C.  dehilis.  and  Sitzschia  sp..  com¬ 
prised  v6  1  t"  „  of  the  cells  examined.  The  remaining  cells 
were  classified  as  centrics.  pennates.  or  dinoflagellates. 
Aggregates  were  very  loosely  held  together  by  tangled 
spines.  Neither  visual  nor  microscopic  examination  revealed 
the  presence  of  abundant  mucus  exudates  within  the 
aggregates.  The  average  cell  volume  of  diatoms  within  the 
floes  w  as  calculated  based  on  the  volumes  of  the  three  major 
diatom  species,  assuming  their  shapes  to  be  cylinders.  A  cell 
with  a  volume  equivalent  to  the  average  volume  of  cells 
within  the  floes  was  calculated  to  have  an  average  radius  of 
7  /im.  This  average  radius  was  used  in  subsequent  calcula¬ 
tions  for  diatom  characteristics. 

The  porosity  of  1 2  diatom  tlocs  is  shown  as  a  function  of 
aggregate  volume  in  Fig.  2.  Aggregate  volumes  ranged  be¬ 
tween  0.19  and  4.2  cm  ':  based  on  these  volumes,  the  diame¬ 
ter  of  a  spherical  aggregate  of  equal  volume  varied  between 
7  and  20  mm.  Porosity  increased  non-linearly  with  aggregate 
volume,  ranging  between  0.99931  and  0.99984.  A  maximum 
error  in  porosity,  .ip.  was  determined  using  a  first-term 


AGGREGATE  VOLUME  cm3 

Fig.  2.  Aggregate  porosity  ip)  for  12  diatom  tlocs  as  a  function  of 
Hoc  volume  (  V.  cm’),  with  regression  line  p  =  1  -3.5  x  I0~4  y~°  i0 


Table  1.  Percentage  composition  of  diatom  floes  collected  from  the 
Santa  Barbara  Channel.  July  1986.  by  phytoplankton  cell  number. 
Reported  percentages  ignore  whether  cells  present  in  chains. 
I.V=  10.  S.D.  =  standard  deviation) 


Species 

Mean  % 

( +  S.D.) 

Mean  cell 
dimensions 
(diameter  x 
height,  /im) 

Mean 

cell 

volume 

(/<m) 

Chaetocerns  debilis 

18  ±2 

21  x  15 

5440 

Chaetoceros  radicans 

58+6 

9x  13 

780 

\itzschia  sp. 

10  +  2 

7x30 

1  150 

Miscellaneous 

Centrics 

8  +  3 

15x15 

Pennates 

3  ±  1 

10x30 

Dinoflagellates 

4  +  2 

12  x  2 

Tay  lor  Series  expansion  (Mickley  et  al.  1957)  for  porosity  as 
a  function  of  three  variables,  or 
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Using  Eq.  1.  this  expression  becomes: 
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To  calculate  J  p.  we  used  measured  values,  and  assumed  the 
following  maximum  errors  based  on  data  in  Alldredge  and 
Gotschalk  (1988):  JtF=+10%.  J  P  =  ±  0.08  g  cm  ’.  and 
11=  ±  30%.  As  shown  in  Table  1.  the  errors  in  porosity  are 
small  and  are  <  0.03°  o. 

Measured  settling  velocities  of  diatom  floes  in  situ 
ranged  between  49  and  190  md  with  an  average  of 
1 10  +  40  m  d  1  (Fig.  3).  We  also  calculated  the  settling  ve¬ 
locities  using  Eqs.  I  to  4.  Predicted  settling  velocities  of  the 
12  diatom  floes  ranged  between  73  and  120  m  d  1  and  were 
not  statistically  different  (Students'  t-test)  from  the  observed 
settling  velocities  of  these  floes,  even  though  the  calculated 
velocities  did  not  appreciably  increase  with  aggregate 
diameter.  The  average  predicted  settling  velocity  of 
100  ±  12  m  d~  1  does  not  differs  significantly  from  the  in  situ 
settling  velocity  of  110  m  d  *.  Similar  values  for  the  ob¬ 
served  and  predicted  settling  velocities  support  the  calcula¬ 
tion  of  an  equivalent  aggregate  radius  used  in  this  analysis, 
and  the  application  of  the  permeable  aggregate  model  for 
diatom  floes. 

Measured  sinking  velocities  were  used  in  Eqs.  5  and  6  to 
calculate  intra-aggregate  flow  velocities  between  24  and 
1 60  pm  s  1  ( 2  to  1 4  m  d  1 ),  with  an  average  intra-aggregate 
velocity  of  70/ims  1  (6md"  ‘).  or  about  5%  of  the  observed 
aggregate  settling  velocity  for  12  diatom  floes  (Table  2).  A 
wide  range  of  intra-aggregate  velocities  of  60  to  160  pm  s" 1 
are  predicted  for  aggregates  in  a  narrow  size  range  of  17  to 
20  mm  in  diameter.  This  reflects  the  variable  porosity,  and 
therefore,  permeability,  of  diatom  floes. 

The  overall  rate  of  nutrient  uptake  by  aggregated  di¬ 
atoms  undergoing  gravitational  settling  was  compared  to: 
(1)  unattached  cells  in  laminar  shear  flow  and  (2)  unattached 
cells  sinking  in  an  undisturbed  fluid,  using  a  Relative  Uptake 
Factor  analysis.  Fluid  shear  rates  within  the  ocean  mixed 
layer  based  on  turbulent  energy  dissipation  rates  are  highly 
variable,  with  ranges  of  0.15  to  0.44s" 1  (Mourn  and  Cald¬ 
well  1985).  0.001  to  Is'1.  (Shay  and  Gregg  1984).  and  0.01 
to  0.04s"1  (Oakey  and  Elliott  1982).  The  larger  the  shear 
rate,  the  greater  the  uptake  of  an  unattached  cell  and  the 
smaller  the  Relative  Uptake  Factor.  In  the  first  case,  we 
assumed  the  largest  fluid  shear  rate  of  1  s" 1  as  an  upper  limit 
of  fluid  shear  for  unattached  cells.  This  results  in  a  conserva¬ 
tive  evaluation  of  the  benefits  of  diatom  nutrient  uptake 
within  a  settling  aggregate.  For  the  12  diatom  floes,  we 
calculated  Relative  Uptake  Factors  between  1.0  and  2.1 
(average  of  1.4),  indicating  that  diatoms  within  the  aggregate 
could  potentially  utilize  dissolved  nutrients  up  to  2.1  times 
faster  than  unattached  diatoms  experiencing  fluid  shear  in 
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Fig.  3.  Observed  and  calculated  settling  velocities  of  diatom  Hoes 


Table  2.  Properties  of  diatom  lloes  and  predicted  intra-aggregate 
velocities  and  relative  uptake  factors 


Hoc 

diam. 

( mm  1 

Hoc  porosity 

Measured 
settling 
velocity 
i  m  d  1 ) 

Intra- 
aggregate 
velocity 
turn  s ' 1 ) 

Relative 

uptake 

factor 

71 

(144434  -  0.01)029 

49 

v- 

1.3 -0.3 

K)  > 

(1. 9943S  -0.00029 

X9 

36 

1.2  -0.4 

Kill 

0  44443  -0.0002' 

62 

24 

1 .0  -  0.4 

111) 

0.44959  -  0  00019 

50 

26 

1.1  -  0.3 

1 2  11 

0.49960  -  0.00019 

46 

45 

13  -0.3 

14.(1 

0.44974  -0.00012 

100 

71 

1.6  -0.3 

15.0 

0.99965  X  0.0001 7 

160 

67 

1.5 +  0.4 

r.o 

0.44974  x  0.00010 

1  10 

70 

1.6  -0.3 

r  ii 

0,999X4  -O.OOOOS 

190 

160 

2.1  +  0.4 

IS, II 

0.994X1  -0.00009 

170 

no 

1.8  -0.4 

20  (1 

0.994X5  x  0.00007 

150 

120 

1.9  -0.4 

20.0 

0.94479  =  0.00010 

130 

60 

1.5  ±0.3 

the  water  column  (Table  2).  This  analysis  accounts  for  nu¬ 
trient  depletion  within  the  aggregate,  and  requires  that  both 
unattached  diatoms  and  the  diatom  floe  are  exposed  to  low 
nutrient  concentrations  where  a  first  order  uptake  model 
accurately  describes  nutrient  uptake  kinetics. 

In  case  2,  Relative  Uptake  Factors  were  also  calculated 
assuming  unattached  cells  sank  in  undisturbed  fluid. 
Smayda  and  Boleyn  (1967)  found  that  the  mean  sinking  rate 
of  chains  of  Chaetoceros  lauderi  up  to  20  cells  long  ranged 
from  0.46  to  1.54  md  with  no  clear  relationship  between 
sinking  speed  and  culture  age.  This  species  had  a  cell  diame¬ 
ter  ranging  from  1 1  to  28 /tm  making  it  similar  in  size  to  the 
two  Chaetoceros  species  that  made  up  an  average  of  76%  of 
the  diatoms  in  our  floes.  Assuming  a  sinking  velocity  of 
1  md  ' 1  for  unattached  cells,  the  Relative  Uptake  Factors  of 
cells  within  aggregates  ranged  between  1.1  and  2.2.  This 
suggests  that  consideration  of  unattached  cells  as  either 
sinking  or  suspended  within  a  laminar  shear  field  is  un¬ 
important  in  obtaining  Relative  Uptake  Factors  greater 
than  unity  for  diatom  floes. 

The  calculation  of  a  Relative  Uptake  Factor  is  based  on 
a  series  of  measured  values  and  equations.  The  maximum 
error  in  this  calculation  can  be  estimated  using  a  Taylor 


series  expansion  about  Uq.  It),  by  calculating  errors  in  o. 
using  Fqs.  1 1  and  12.  and  Shu,  using  Eq.  X.  An  error  for  Sh,, 
was  not  determined  since  an  upper  limit  of  G  was  used  in 
the  calculation  of  this  Sherwood  number.  The  error  in  o 
was  determined  ."or  constant  EH.  D.  a  .  and  t  .  with  the  fol¬ 
lowing  assumed  errors:  !U'=  *  ! 0‘* J P  =  +0.08 gem  \ 
J  I ,  =  ±  30"  n.  and  Jit  =  +  50%.  These  maximum  errors  are 
shown  in  Table  2  for  each  Relative  Uptake  Factor.  On  the 
average.  ;  =  1.4  ±0.3.  which  is  larger  than  unity.  Within 
these  large  error  bounds,  four  Relative  Uptake  Factors  were 
less  than  unity,  indicating  that  within  some  aggiegates.  all 
cells  may  not  have  increased  nutrient  uptake  compared  to 
unattached  cells. 


Discussion 

The  porosities  of  diatom  Hoes  were  highly  variable,  but  were 
in  excess  oft). 999 10  for  the  12  diatom  floes  analyzed.  These 
values  are  higher  than  average  porosities  estimated  for 
aggregates  found  in  engineered  bioreactors.  For  example,  po¬ 
rosities  for  activated  sludge  floes  range  from  0.79  to  0.90. 
(Mueller  et  al.  1966).  and  porosities  for  mold  pellets  range 
from  0.800  to  0.988  (Yano  et  al.  1961).  Since  the  size  of  an 
aggregate  is  a  function  of  fluid  shear  (Hunt  1986).  the  high 
diatom  floe  porosities  are  reasonable  for  aggregates  formed 
in  low  shear  environments  typical  of  natural  systems.  In 
engineered  systems,  such  as  wastewater  treatment  systems, 
shear  rates  are  in  the  range  of  100  to  1000s" l.  Uarge  shear 
rates  contribute  to  aggregate  breakup,  re-compaction,  and 
the  formation  of  smaller  aggregates  with  reduced  porosities. 
Porosities  of  diatom  lloes  are  similar  to  aggregate  porosities 
determined  using  collision  calculations  based  on  random 
particle  trajectories  in  computer  models  (Tambo  and  Wata- 
nabe  1979).  These  computer-generated  aggregates,  de¬ 
scribed  as  fractals  (Witten  anc  Cates  1986).  typically  neglect 
aggregate  breakup. 

The  Relative  Uptake  Factors  for  the  12  diatom  floes 
were  determined  from  several  aggregate  characteristics,  and 
a  series  of  relationships  drawn  from  the  literature,  to  be  in 
the  range  of  1.0  +  0.3  to  2.1  ±  0.4  (average  1.4  +  0.3).  This 
indicates  flocculation  can  increase  nutrient  uptake  by  at¬ 
tached  diatoms.  In  our  error  analysis,  we  estimated  the 
cumulative  effect  of  errors  in  measured  parameters  still  re¬ 
sulted  in  Relative  Uptake  Factors  greater  than  unity.  Two 
equations  were  not  included  in  the  error  analysis:  Eq.  5.  used 
to  relate  aggregate  porosity  and  permeability,  and  Eq.6. 
usrd  to  determine  the  intra-aggregate  velocity.  Several  cor¬ 
relations  have  been  found  to  provide  good  agreement  be¬ 
tween  experimentally  determined  permeabilities  and  porosi¬ 
ties  for  fibrous  media  at  Re  <  10  (Jackson  and  James  1986). 
However,  we  are  unaware  of  correlations  for  high-porosity 
aggregates  at  the  slightly  larger  Reynolds  numbers 
(3  <  Re  <  31 )  determined  for  settling  diatom  floes.  As  pre¬ 
viously  stated,  we  used  the  Davies  correlation  since  this 
expression  was  supported  by  experiments  conducted  by 
Masilyah  and  Polikar  (1980)  for  0.2  <  Re  <  120.  Instead  of 
performing  an  error  analysis  on  Eqs.  5  and  6.  we  calculated 
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the  error  in  the  Relative  Uptake  Factor  by  assuming  the 
intra-aggregate  velocity  obtained  from  these  equations  was 
only  accurate  within  ±  5 0“u.  Even  with  this  large  variation 
in  intra-aggregate  velocity,  only  four  of  the  Relative  Uptake 
Factors,  within  a  standard  error,  were  less  than  unity. 
Therefore,  while  our  hypothesis  remains  physically  un¬ 
tested.  our  calculations  provide  strong  evidence  that  diatom 
flocculation  can  increase  nutrient  uptake. 

There  are  two  defined  limits  to  a  Relative  Uptake  Factor 
analysis  of  mass  transfer  to  attached  diatoms.  The  lower 
limit  results  if  the  advective  How  of  nutrients  into  the 
aggregate  is  less  than  the  diffusive  flux.  This  occurs  when  the 
Peclet  number  for  the  Hoc.  Pe;  =  </,  u  D  <<  1  (Logan  and 
Hunt  1988).  For  the  diatom  Hoes  in  this  study,  this  would 
occur  for  it  <?  0.01  md  ‘.a  velocity  much  lower  than  calcu¬ 
lated  values  (  Table  2).  However,  since  unattached  cells  settle 
at  '■  I  m  d  *.  intra-aggregate  flows  would  have  to  be  greater 
than  1  md  1  to  confer  any  additional  nutrient  advantage  to 
attached  diatoms. 

The  upper  limit  of  the  Relative  Uptake  factor  is  defined 
by  the  ratio  of  Sherwood  numbers  in  Eq.  10.  As  the  intra¬ 
aggregate  velocity  increases,  and  the  aggregate  porosity  de¬ 
creases.  the  ratio  (1  —  e'*]  <p  approaches  unity.  Therefore, 
the  overall  uptake  by  attached  cells  cannot  exceed  the  up¬ 
take  by  a  cell  fixed  in  a  flow  field  of  velocity  u.  divided  by 
the  uptake  of  a  cell  in  laminar  shear  (or  sinking  at  1  md  ' l). 
However,  this  ratio  of  Sherwood  numbers  may  still  under¬ 
estimate  uptake  by  attached  cells  since  the  high  fluid  veloci¬ 
ties  at  an  aggregate  surface  are  not  included  in  our  calcula¬ 
tions.  We  assumed  the  Hoc  was  homogeneous  and  spherical 
in  shape,  and  ttiat  all  cells  within  the  aggregate  resided  in  a 
fluid  field  described  by  the  average  intra-aggregate  fluid 
velocity.  This  oversimplification  of  the  structure  and  shape 
of  a  marine  snow  aggregate  was  necessary  to  obtain  an 
analytical  solution.  Natural  aggregates  often  have  long 
comet-shaped  tails  and  irregular  shapes:  the  exact  shape  of 
the  aggregate,  however,  is  not  expected  to  be  important  in 
mass  transfer  calculations  (Aris  1975).  More  important  to 
our  calculations  is  the  assumption  of  the  homogeneity  of  the 
aggregate.  Many  aggregates  contain  internal  voids  (see  All¬ 
dredge  and  Silver  1988,  Fig.  2a,  and  Alldredge  and  Got- 
schalk  1988,  Fig.  1  b  for  visual  examples).  As  a  result,  many 
cells  (more  than  for  a  sphere)  will  be  at  the  aggregate  surface 
and  will  experience  fluid  velocities  much  larger  than  the 
average  intra-aggregate  velocity  calculated  for  the  homoge¬ 
neous  spherical  aggregates  in  this  study.  These  surface  cells 
could  remove  nutrients  much  faster  than  cells  within  the 
aggregate  interior.  This  effect  would  increase  the  benefits  of 
attachment. 

The  maximum  benefit  of  attachment  for  diatoms  in  per¬ 
meable  aggregates,  for  any  model,  is  limited  by  the  maxi¬ 
mum  possible  increase  in  uptake  with  respect  to  fluid  hydro¬ 
dynamics.  This  hydrodynamic  effect  is  described  mathemat¬ 
ically  in  terms  of  Sherwood  number  correlations  for  cells  in 
laminar  fluid  shear  and  for  spheres  fixed  in  a  fluid  field. 
Despite  work  by  Munk  and  Riley  (1952)  over  three  decades 
ago.  little  improvement  has  been  made  in  quantifying  the 
effect  of  hydrodynamics  on  nutrient  uptake  by  phytoplank¬ 
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ton.  The  physics  of  mass  transfer  to  spheres  and  other  ob¬ 
jects  is  well  understood  in  the  field  of  chemical  engineering, 
and  kinetic  models  relating  enzyme  kinetics  to  bulk  nutrient 
concentration  are  well  studied  by  biologists.  However,  the 
increase  in  uptake  kinetics  as  a  result  of  fluid  hydrodynam¬ 
ics  in  biological  systems  is  not  well  understood. 

To  our  knowledge,  only  two  studies  on  the  effect  of 
hydrodynamics  on  nutrient  uptake  by  marine  diatoms  can 
be  used  to  calculate  a  Sherwood  number.  The  first  study,  by 
Canelli  and  Fuhs  ( 1976).  investigated  phosphorus  uptake  by 
cells  of  Thalassiosini  Jluviatilis  fixed  in  a  flow  field  and  is 
used  in  our  calculations  (Eq.8).  In  Fig. 4.  two  other  Sher¬ 
wood  number  correlations  that  have  been  used  for  similar 
purposes  in  other  studies  are  compared  to  Eq.8.  The  first 
correlation,  developed  by  Brian  and  Hales  (Sherwood  et  al. 
1975).  and  used  by  Logan  and  Hunt  (1987.  1988)  to  describe 
uptake  by  cells  fixed  in  a  flow  field  is: 

Shu  =  <1  +  0.48Pe;  V  ’.  (15) 

A  second  correlation  was  used  by  Munk  and  Riley  (1952)  in 
their  classic  study  for  sinking  spheres;  their  correlation  is: 

Sh  =  1  +  0.5  Pe  4-  0.60  Pe2 .  (16) 

Fig. 4  shows  that  at  velocities  <160 pirns'1  both  Eqs.  15 
and  16  underestimate  the  effect  of  fluid  flow  when  compared 
to  the  Sherwood  number  correlation  derived  using  the  ex¬ 
perimental  data  of  Canelli  and  Fuhs.  All  three  correlations 
indicate  Sherwood  numbers  should  continue  to  increase 
with  fluid  velocities  greater  than  160/uns~ ',  but  this  obser¬ 
vation  is  not  supported  by  Canelli  and  Fuhs'  data  (Fig.  1). 
These  results  suggest  there  is  some  effect  of  fluid  velocity 
that  is  not  adequately  described  in  our  mass  transfer  analy¬ 
sis.  Since  the  calculated  intra-aggregate  velocities  shown  in 
Table  2  are  less  than  <  160 /tms'1,  the  Relative  Uptake 
Factors  in  Table  2  are  conservative. 

The  second  study  that  can  be  used  to  determine  a  Sher¬ 
wood  number  correlation  was  reported  by  Pasciak  and 
Gavis  (1975,  Fig.  5)  on  nitrate  uptake  by  Ditylum  brightwel- 
li.  Shown  in  Fig.  5  is  a  comparison  of  Sherwood  numbers 
for  a  41  ft m  cell  (the  equivalent  spherical  radius  of  D.  bright- 
wetli )  to  Sherwood  numbers  calculated  using  Eq.9.  Much 
greater  uptake  for  unattached  diatoms  as  a  result  of  fluid 
shear  is  predicted  from  Eq.  9  than  was  observed  by  Pasciak 
and  Gavis  (1975).  This  suggests  that  we  overestimated  the 
uptake  of  unattached  diatoms  in  our  calculations,  which 
would  lead  to  a  larger  Relative  Uptake  ratio  than  shown  in 
Table  2.  Thus,  experimental  work  by  Pasciak  and  Gavis 
(1975),  as  well  as  by  Canelli  and  Fuhs  (1976),  suggests  that 
our  comparison  of  the  nutrient  advantage  experienced  by 
aggregated  over  unattached  diatoms  is  conservative. 

Flocculation  has  historically  been  viewed  as  a  disadvant¬ 
age  to  attached  cells  for  two  reasons.  First,  aggregates  have 
been  viewed  as  impermeable,  resulting  in  decreased  nutrient 
availability  to  cells  within  the  aggregate  interior  (Yano  et  al. 
1961).  Second,  since  flocculation  increases  settling  velocity, 
diatoms  may  be  removed  faster  from  the  euphotic  zone  than 
unattached  cells.  Smetacek  (1985),  for  example,  has  argued 
that  aggregation  and  mass  sinking  represent  a  transition 
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FLUID  VELOCITY  gm  s'" 


Fig.  4.  Sherwood  numbers  of  a  eell  7  /< m  in  radius  Fixed  in  a  flow 
Field,  ealeulated  using:  F.q.  8.  based  on  data  in  Canelli  and  Fuhs 
t  I97ft)(dashed  line);  Eq.  16.  the  Brian  and  Hales  Correlation  (Sher¬ 
wood  et  al.  1975)  (eontinuous  line):  Eq.  15.  a  correlation  used  by 
Munk  and  Riley  (1957)  (doited  line) 
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Fig.  5.  Sherwood  numbers  calculated  using  Eq.  9  (dashed  line)  for 
a  cell  41  fim  in  radius  as  a  function  of  fluid  shear  rate  versus  results 
adapted  from  Pasciak  and  Gavis  (1975)  (continuous  line) 


from  a  growing  to  a  resting  stage  in  the  life  histories  of 
diatoms.  Aggregated  cells  sink  out  of  nutrient-depleted  sur¬ 
face  waters  into  deep,  nutrient-rich  water  and  await  re- 
suspension  back  into  the  euphotic  zone  as  resting  stages. 
While  this  evolutionary  explanation  for  the  mass  floccula¬ 
tion  of  diatom  blooms  may  be  correct,  resting  stages  are  rare 
in  diatom  floes  observed  in  surface  waters  (Alldredge  and 
Gotschalk  in  press)  although  they  may  form  at  depth.  Our 
data  suggest  that  the  flocculation  of  diatom  blooms  may  be 
immediately  adaptive  because  aggregation  alters  the  fluid 
environment  of  a  cell.  A  single  unattached  cell  or  chain  is 
contained  within  a  microscale  eddy  and  must  move  with  the 
bulk  fluid;  however,  fluid  flow  around  and  through  a  rapid¬ 
ly  sinking  aggregate  actually  alters  the  fluid  environment  of 
associated  cells,  and  increases  nutrient  uptake  by  attached 
cells  compared  to  unattached  cells.  Flocculation  may,  there¬ 
fore,  be  advantageous  to  nutrient-stressed  diatoms  by  in¬ 
creasing  the  potential  for  uptake  of  scarce  nutrients  follow¬ 
ing  an  intense  phytoplankton  bloom. 

Evolutionary  theory  predicts  that,  ail  else  being  equal, 
the  genotype  producing  aggregation  behavior  will  be  se¬ 
lected  for  only  if  diatoms  carrying  that  genotype  remain  in 


the  mixed  layer  long  enough  to  experience  differential 
growth  rates  above  those  of  non-aggregating  genotypes  of 
the  same  species.  If  the  magnitude  of  this  differential  repro¬ 
duction  is  very  small,  more  time  will  be  required  to  establish 
the  genotype  for  aggregation  within  the  species.  However,  as 
long  as  the  differential  reproduction  of  aggregating  diatoms 
is  greater  than  non-aggregating  cells,  i.e..  the  fitness  of  the 
aggregation  genotype  is  greater  than  unity,  then  aggregation 
behavior  will  be  selected  for  and  the  aggregation  genotype 
will  come  to  dominate  the  population  (see  Emlen  1973  for 
discussion).  Phytoplankton  cell  growth  and  reproduction  is 
directly  related  to  nutrient  uptake.  Our  finding  that  aggre¬ 
gated  cells  may  take  up  nutrients  as  much  as  2.1  times  faster 
than  unaggregated  cells  of  the  same  species  suggests  that  the 
differential  production  of  aggregated  cells  has  the  potential 
to  be  quite  high.  However,  even  Relative  Uptake  Factors 
only  slightly  greater  than  unity  could  result  in  gradual  selec¬ 
tion  for  diatoms  which  aggregate. 

We  propose  that  some  aggregating  diatoms  remain  in  the 
mixed  layer  long  enough  to  experience  differential  reproduc¬ 
tion  because  of  two  sequential  processes.  First,  despite  their 
rapid  sinking  rates,  some  diatom  floes  remain  suspended  in 
the  mixed  layer  for  many  days  before  they  sink  out  of  it. 
Lande  and  Wood  (1987)  calculated  that  slowly  sinking  par¬ 
ticles  may  make  many  vertical  excursions  within  the  mixed 
layer  before  they  are  lost  via  settlement.  Empirical  evidence 
indicates  that  some  rapidly  settling  particles  also  make  such 
excursions.  For  example.  5  to  10%  of  macro-crustacean 
fecal  pellets  sinking  at  rates  of  18  to  170md*‘  have  been 
observed  suspended  in  the  mixed  layer  10  d  after  produc¬ 
tion  (Alldredge  et  al.  1987).  Moreover,  diatom  floes  several 
days  old  do  indeed  occur  abundantly  in  surface  waters  off 
southern  California  (Alldredge  and  Gotschalk  in  press).  Our 
calculations  indicate  that  cells  attached  to  these  floes,  due  to 
increased  nutrient  uptake,  should  have  a  selective  advantage 
over  nutrient  depleted,  unattached  cells  remaining  an  equal 
amount  of  time  in  the  mixed  layer.  Second,  once  replete,  at 
least  some  cells  must  fall  off  floes  before  the  floes  leave  the 
mixed  layer.  Diatom  floes  are  fragile  and  erode  or  fragment 
at  shear  rates  normally  encountered  in  the  ocean  (McGilli- 
vary  and  Alldredge  1987).  Single  chains  or  floe  fragments 
will  more  easily  become  suspended  than  the  larger  and  heav¬ 
ier  parent  floes.  A  cycle  of  repeated  resuspension,  coupled 
with  gradual  erosion  or  fragmentation  of  only  a  small  per¬ 
centage  of  diatom  floes  formed  after  a  bloom,  would  be 
adequate  to  reseed  the  mixed  layer  with  unattached  vegeta¬ 
tive  cells  replete  with  nutrients  due  to  their  prior  association 
with  the  floes.  It  is  these  cells  which  may  be  capable  of 
sustained  survival  following  a  bloom  and  which  eventually 
re-establish  diatom  populations  in  surface  waters  when 
more  favorable  nutrient  conditions  reoccur. 
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Fractal  geometry  of  marine  snow  and  other 
biological  aggregates 


Abstract — Fractal  dimensions  of  aggregates  can 
potentially  be  used  to  classify  aggregate  mor¬ 
phology  as  well  as  to  identify  coagulation  mech¬ 
anisms.  Microbial  aggregates  of  Zoogloea  ramig - 
era  have  a  cluster  fractal  dimension  of  1.8 ±0.3 
(±SD),  suggesting  that  these  aggregates  are  formed 
through  cluster-cluster  coagulation.  An  analysis 
of  size-porosity  correlations  for  two  types  of  ma¬ 
rine  snow  aggregates  yielded  fractal  dimensions 
of  1.39±0.06  and  1.52±0.19,  which  were  lower 
than  values  describing  inorganic  colloidal  aggre¬ 
gation. 


Large,  amorphous  aggregates,  known  as 
marine  snow,  compose  the  largest  compo¬ 
nent  of  mass  sedimenting  through  the  water 
column  (Fowler  and  Knauer  1986).  These 
aggregates  consist  of  diatoms,  bacteria,  fecal 
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pellets,  cast  houses  ofappendicularians,  and 
nearly  all  other  microscopic  organic  and  in¬ 
organic  matter  present  in  the  ocean.  The 
manner  in  which  these  larger  aggregates  are 
formed  is  not  well  understood.  The  growth 
of  microaggregates  (<0.5  mm)  and  marine 
snow-sized  aggregates  (>0.5  mm)  can  result 
from  Brownian  motion,  eddy  diffusion, 
shear  coagulation,  and  differential  sedimen¬ 
tation,  as  well  as  through  growth  of  organ¬ 
isms  within  the  aggregate  microhabitat 
McCave  (1984)  has  calculated  that  shear 
rates  in  the  ocean  are  too  low  to  account  for 
observed,  steady  state  particle  size  distri¬ 
butions.  Although  this  conclusion  implies 
that  physical  mechanisms  of  aggregate  for¬ 
mation  are  less  important  than  biological 
mechanisms  in  marine  systems,  marine 
snow  aggregates  have  macroscopic  mor¬ 
phologies  that  are  characteristic  of  inorganic 
aggregates  formed  through  physical  coagu¬ 
lation  mechanisms.  This  similarity,  which 
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can  be  established  with  fractal  geometry 
(Mandelbrot  1977;  Feder  1988),  provides 
evidence  that  aggregates  can  form  in  the 
ocean  through  physical  coagulation. 

The  structures  of  inorganic  aggregates, 
such  as  gold  panicle  dusters  formed  irre¬ 
versibly  in  suspension,  are  scale  invariant 
or  fractal  (Witten  and  Cates  1986).  The 
number  of  panicles,  N,  in  a  fractal  aggregate 
is 

N  -  /°"  (1) 

where  D„  is  the  duster  fractal  dimension 
determined  for  the  object  in  n  dimensions, 
and  /  is  the  characteristic  length  scale  of  the 
aggregate.  For  a  Euclidean  object  such  as  a 
sphere,  Z)3  =  3.  Aggregates  formed  by  ran¬ 
dom  processes  have  fractal  dimensions  sig¬ 
nificantly  less  than  the  Euclidean  3.  Com¬ 
puter  simulation  of  aggregate  growth  by 
random  processes  indicates  that  the  mag¬ 
nitude  of  the  fractal  dimension  is  deter¬ 
mined  by  the  mechanism  of  aggregate 
growth.  Aggregates  formed  through  the  ad¬ 
dition  of  particles  into  the  cluster  one  at  a 
time  (particle-cluster)  have  three-dimen¬ 
sional  fractal  dimensions  in  the  range  of  2.5- 
3.0  (Schaefer  1989).  Aggregates  formed 
through  collision  of  dusters  (cluster-cluster) 
have  lower  fractal  dimensions,  typically  with 
1.6  <  D3  <  2.2  (Witten  and  Cates  1986) 

Recent  research  on  inorganic  colloidal  ag¬ 
gregation  has  demonstrated  a  relationship 
between  aggregate  structure  and  cluster- 
cluster  aggregation  kinetics  (Lin  et  al.  1989). 
Computer  simulations  and  experiments  with 
several  types  of  colloids  indicate  that  two 
different  types  of  aggregates  form  as  a  func¬ 
tion  of  different  particle  stickiness  or  dif¬ 
ferent  collision  efficiencies  (Table  1).  As  the 
probability  of  attachment  on  collision  ap¬ 
proaches  unity,  particles  stick  upon  contact, 

■  forming  highly  tenuous  structures.  This  pat¬ 
tern  results  in  very  rapid  coagulation,  re¬ 
ferred  to  as  diffusion-limited  colloidal  ag¬ 
gregation  (DLCA)  because  the  rate  of 
aggregate  growth  is  limited  by  particle  trans¬ 
port  The  highly  tenuous  structure  is  a  result 
of  very  sticky  particles  moving  in  a  random 
walk;  particles  cannot  penetrate  the  aggre¬ 
gate  since  they  attach  to  particles  with  which 
they  collide  on  the  aggregate  exterior.  Ag¬ 
gregates  formed  through  DLCA  are  char- 


Table  1.  Fractal  dimensions  of  inorganic  colloids 
formed  through  either  DLCA  or  RLCA  (data  from  Lin 
et  al.  1989). 


Colloid 

DLCA 

RLCA 

Gold 

1.86 

2.14 

Silica 

1.85 

2.07 

Polystyrene 

1.82 

2.09 

acterized  by  fractal  dimensions  around  1.8. 
As  the  attachment  probability  approaches 
zero,  particles  may  collide  many  times  be¬ 
fore  sticking,  allowing  them  to  penetrate  the 
aggregate  and  increasing  the  aggregate  den¬ 
sity.  This  aggregation  is  classified  as  reac¬ 
tion-limited  colloidal  aggregation  (RLCA) 
and  results  in  aggregates  with  D3  values 
around  2.1. 

The  occurrence  of  two  distinct  regimes  of 
colloidal  aggregation  suggests  that  fractal 
scaling  relationships  can  potentially  identify 
the  aggregate  formation  mechanism  as  well 
as  the  attachment  efficiency  of  particles 
composing  the  aggregate.  Lin  et  al.  (1989) 
demonstrated  universality  of  fractal  classi¬ 
fication  for  inorganic  systems  with  colloidal 
gold,  silica,  and  polystyrene.  It  has  not  been 
established,  however,  that  biological  aggre¬ 
gates  can  be  similarly  classified.  Most  mi¬ 
croorganisms  grown  in  batch  culture  aggre¬ 
gate  during  late  log  growth  (Calleja  1984). 
This  aggregation  creates  a  more  complex 
system  than  previously  examined  for  inor¬ 
ganic  coagulation  because  particles  (daugh¬ 
ter  cells)  can  be  produced  during  coagula¬ 
tion.  Microbial  growth  occurs  during 
aggregation,  so  we  would  expea  biological 
aggregates  to  have  D3  values  that  are  be¬ 
tween  values  reported  for  particle-cluster 
and  cluster-cluster  models.  The  size  of  the 
microbial  aggregates  and  the  range  of  co¬ 
agulation  mechanisms  (e.g.  shear,  differen¬ 
tial  sedimentation),  however,  create  con¬ 
ditions  far  different  from  systems  previously 
characterized  by  fractal  analysis. 

In  order  to  establish  that  biological  ag¬ 
gregates  have  fractal  structures  similar  to 
inorganic  aggregates,  we  developed  micro¬ 
bial  aggregates  composed  of  homogeneous¬ 
ly  sized  particles  with  pure  cultures  of  mi¬ 
croorganisms.  We  grew  aggregates  of 
Zoogloea  ramigera  type  I-16-M  (ATCC- 
19623),  a  rod-shaped,  gram-negative  bac- 
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Fig.  1 .  Cells  per  floe  of  Zoogloea  ramigera  (■)  with 
a  regression  line  through  the  data  (r2  =  0.87)  indicating 
D j  =  1.8 ±0.3  (±  1  SD).  Data  points  represent  averages 
of  3-6  aggregates. 


terium  which  attaches  with  cellulosic  fibers 
that  resist  floe  breakage  during  aggregate- 
sizing  experiments.  Suspended  cultures  (15 
ml)  were  developed  at  room  temperature 
(23°C)  on  nutrient  broth  (Difco)  in  36-ml 
test  tubes  clamped  to  a  laboratory  rotator. 
After  24  h,  cultures  were  poured  into  a  ster¬ 
ile  Petri  dish  and  diluted  with  sterile  water 
to  increase  floe  separation.  Individual  ag¬ 
gregates  of  Z.  ramigera  were  captured  in 
~20  n\  of  media  with  a  100-m1  pipette  with 
the  plastic  tip  cut  between  the  25-  and  50- 
/zl  markings  and  transferred  to  a  2-ml  drop 
of  solution  to  aid  separation  of  floes  from 
free-living  cells.  The  pipet  tip  was  then  rinsed 
with  sterile  water  and  the  aggregate  trans¬ 
ferred  to  a  well  slide.  Floes  were  viewed  with 
direct  light  microscopy  (10  x  power,  Olym¬ 
pus  model  BH-2)  to  determine  the  largest 
aggregate  length  and  to  sketch  the  aggregate. 

The  number  of  cells  composing  each  mi¬ 
croscopically  sized  aggregate  was  deter¬ 
mined  by  rinsing  the  slide  contents  into  a 
sterile,  acid-washed  vial  and  adding  sterile 
water  for  a  final  volume  of  1  ml.  The  floe 
was  dispersed  enzymatically  with  10  mg  ml-1 
of  cellulase  (Sigma  Chemical  No.  C-7377), 
vortexing  the  solution  for  1  min,  and  plac¬ 
ing  the  vials  in  a  rotator.  Samples  were  pre¬ 
served  with  2%  formaldehyde  and  counted 
with  acridine  orange  epifluorescence  (Hob- 
bie  et  al.  1977).  Samples  containing  incom¬ 
pletely  dispersed  floes  were  discarded. 


Fig.  2.  Occupied  volume  (1  -  porosity)  correla¬ 
tions  with  aggregate  size  for  several  different  types  of 
biological  aggregates  from  both  natural  and  engineered 

environments:  Zoogloea  ramigera  ( - );  marine  snow 

( - );  diatom  floes  ( — );  normal  wastewater  bio¬ 

reactor  floes  ( );  filamentous  wastewater  bioreactor 

floes  ( - ). 

Transfer  of  culture  material  without  floes 
through  the  above  sequence  of  steps  was 
used  to  verify  that  free-living  cells  included 
with  the  floe  transfer  were  a  negligible  por¬ 
tion  of  the  total  counts. 

Fractal  aggregates  are  characterized  by 
large  spatial  and  sample-to-sample  fluctua¬ 
tions  in  properties  (Meakin  1988).  There¬ 
fore,  we  sized  47  aggregates  by  longest  length 
and  averaged  all  aggregates  within  0.1 -mm 
size  classes  spanning  0.7-1 .5  mm  long.  From 
a  log-log  regression  with  Eq.  1  and  data 
shown  in  Fig.  1,  we  determined  a  fractal 
dimension  for  Z.  ramigera  floes  of  1.8  with 
a  standard  error  of  0.3  (r2  *  0.87,  n  *  8) 
based  on  the  longest  aggregate  length.  This 
value  is  much  smaller  than  the  Euclidean 
value  of  D3  «  3  expected  for  spherical  ag¬ 
gregates  and  is  experimental  evidence  of  the 
three-dimensional  fractal  structure  of  these 
aggregates.  The  magnitude  of  this  fractal  di¬ 
mension  is  within  the  range  indicated  for 
aggregates  formed  through  cluster-cluster 
coagulation,  but  is  not  sufficiently  precise  to 
identify  whether  cluster  formation  is  within 
the  range  expected  for  diffusion-  or  reac¬ 
tion-limited  aggregation. 

Several  investigators  have  related  the  po- 
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Table  2.  Fractal  dimensions  of  biological  aggregates  determined  from  empirical  size-porosity  equations. 


Biol. 

segregate  type 

Sue  range 
(mm) 

Fractal 

dimension 

Reference* 

Marine  snow 

General 

0.4-20 

1.39+0.15* 

Alldredge  and  Gotschalk  1988 

Diatom 

7-20 

1.52±0.19* 

Logan  and  Alldredge  1 989 

Bioreactor 

Normal 

0.2-2 

1.3 

Tambo  and  Watanabe  1979 

Filamentous 

0.3-2 

1.0 

Tambo  and  Watanabe  1979 

Zoogloea  ramigera 

0.4-1. 4 

1. 8  +  0.3* 

This  study 

•*isd 


rosity,  p,  of  biological  aggregates  to  aggre¬ 
gate  size  (Tambo  and  Watanabe  1979;  All¬ 
dredge  and  Gotschalk  1988;  Logan  and 
Alldredge  1 989)  with 

1  -  p  =  alb  (2) 

where  a  and  b  are  empirical  constants,  and 
b  is  negative.  Porosity  has  been  determined 
either  directly  from  gravimetry,  assuming 
an  average  particle  density  (Alldredge  and 
Gotschalk  1988;  Logan  and  Alldredge  1989) 
or  indirectly  from  settling  experiments 
(Tambo  and  Watanabe  1979).  If  the  num¬ 
ber  of  particles  or  cells  comprising  the  ag¬ 
gregate  is  known,  the  porosity  can  be  di¬ 
rectly  calculated  from 


where  Vc  is  the  average  volume  of  cells  com¬ 
posing  the  aggregate  and  V  the  volume  oc¬ 
cupied  by  an  aggregate.  Following  conven¬ 
tion  (Witten  and  Cates  1986;  Li  and 
Ganczarczyk  1987;  Logan  and  All  dredge 
1989),  we  defined  the  volume  of  Z.  ramig- 
era  aggregates  as  the  volume  of  a  sphere  just 
capable  of  enclosing  the  aggregate.  Combin¬ 
ing  Eq.  1  and  3,  we  obtained  the  fractal 
relationship 

(1  -  P)  ~  (4) 

A  comparison  of  Eq.  2  and  4  shows  that  D3 
*3  +  6.  Therefore,  the  fractal  dimension 
of  marine  snow,  as  well  as  other  types  of 
biological  aggregates,  can  be  estimated  with 
published  size-porosity  relationships. 

Although  marine  snow  aggregates  can  be 
much  larger  than  other  types  of  biological 
aggregates,  they  are  denser  than  aggregates 
of  similar  size  (Fig.  2).  The  porosities  de¬ 


termined  for  Z.  ramigera  aggregates  by  our 
microscopic  methods  are  similar  to  poros¬ 
ities  determined  for  aggregates  normally 
found  in  wastewater  treatment  bioreactors. 
On  occasion,  filamentous  cells  become  nu¬ 
merically  abundant  in  these  bioreactors. 
This  produces  highly  filamentous  floes  that 
are  morphologically  different  from  normal 
bioreactor  floes  that  consist  primarily  of 
clumps  of  spherical  or  rod-shaped  cells. 
These  filamentous  floes  are  much  more  po¬ 
rous  than  normal  bioreactor  floes. 

The  fractal  dimensions  for  both  naturally 
occurring  and  engineered  reactor  aggregates 
are  in  the  range  of  1.0-1. 5  (Table  2).  Zoo- 
gloea  ramigera  aggregates  grown  in  the  lab¬ 
oratory  had  the  largest  fractal  dimension  of 
all  aggregates  examined.  This  result  indi¬ 
cates  that  aggregates  grown  in  the  laboratory 
are  more  compact,  having  higher  fractal  di¬ 
mensions  than  aggregates  formed  in  the 
ocean. 

Fractal  dimensions  have  been  calculated 
with  settling  data  (Li  and  Ganczarczyk 
1989).  These  fractal  dimensions,  however, 
contain  a  mixture  of  scaling  properties.  From 
a  force  balance  on  a  settling  particle,  we 
have  (Bird  et  al.  1960) 

V(P  ~  P/)g  =  0.5APjCdIP  (5) 

where  p  is  the  floe  density,  p/the  fluid  den¬ 
sity,  g  the  gravitational  constant,  A  the  pro¬ 
jected  surface  area  of  the  floe  perpendicular 
to  direction  of  settling,  CD  a  drag  coefficient, 
and  U  the  aggregate  settling  velocity.  For 
small  aggregates,  the  drag  coefficient  is  pro¬ 
portional  to  aggregate  settling  velocity  and 
inversely  proportional  to  aggregate  length. 
For  a  solid  object,  such  as  a  sphere  with 
diameter  d,  CD  *  24/Re  for  Re  *  Ud/v  c 
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Table  3.  Fractal  dimensions  of  biological  aggregates  determined  from  settling  velocity  data. 


Biol.  aggregate  type 

Size  range 
(mm) 

Fractal  dimension 

References 

Mannc  snow  (in  situ) 

7-20 

1.26+0.06* 

Alldredge  and 
Gotschalk  1988 

Estuarine 

0.02-2 

1.78 

Gibbs  1985 

Lacustnne 

0.012-0.04 

1.39-1.69 

Hawley  1982 

Oceanic 

-0.1 

2.14 

Oceanic 

~0.4 

1.94 

Recoaguiated  oceanic  sediments 

0.1-1 

1.57 

Kajihara  1971 

•±SD. 


1  where  v  is  the  fluid  kinematic  viscosity. 
Using  this  drag  coefficient  in  Eq.  5,  we  ob¬ 
tain 


V(P  ~  Pr)d 
llApfV 


(6) 


(p  -  pf)  =  (1  -  p)(pc  -  Pf) 

where  pc  is  the  density  of  the  particles  com¬ 
prising  the  aggregate,  so  we  can  define  a 
proportionality  between  settling  velocity  and 
aggregate  geometry  as 


U 


VQ  -  P)d 
A 


(7) 


For  a  spherical,  solid  object  with  Re  1, 
U  ~  d2,  in  agreement  with  Stokes’  law  (Bird 
et  al.  1960). 

We  can  define  a  proportionality  for  a  frac¬ 
tal  object  by  combining  Eq.  1,  3,  and  7,  to 
obtain 


U  ~  p3+l-i>2  (g) 

where  we  have  assumed  that  A  ~  l02,  with 
D2  the  two-dimensional  fractal  dimension 
that  relates  the  projected  surface  area  to  the 
largest  aggregate  length.  Li  and  Ganczar- 
czyk  (1989)  have  defined  the  scaling  power 
in  Eq.  8  as  a  single  fractal  dimension.  This 
scaling  power  is  a  function  of  aggregate  ge¬ 
ometry  in  both  two  and  thee  dimensions 
and  is  defined  here  as  a  fractal  dimension, 
Z>( 2),  without  a  subscript  The  number  in 
parentheses  indicates  the  value  of  this  frac¬ 
tal  dimension  for  a  solid  object  The  fractal 
dimension  D3  can  be  estimated  with  settling 
velocity  data  from  D3  -  D( 2)  -  1  +  D2, 
and  by  assuming  that  D2  =*  2. 

Values  of  D3  obtained  from  Eq.  8  and 
settling  data  span  a  wide  range  of  1 .26-2. 1 4 
(Table  3).  Many  of  these  values  are  below 


those  expected  for  cluster-cluster  coagula¬ 
tion.  Our  estimate  of  Z)3  is  exact  for  any 
self-similar  object  (Euclidean  or  fractal)  if  a 
geometric  length  scale  is  used  to  define  the 
aggregate  area.  If  longest  length  is  used,  then 
D2  ^  2,  and  the  calculated  value  of  D3  is 
an  upper  estimate  of  the  true  value.  Other 
factors  may  influence  the  relationship  be¬ 
tween  D3  and  D( 2).  For  example,  the  marine 
snow  aggregates  examined  by  Alldredge  and 
Gotschalk  (1988)  extend  beyond  the  creep¬ 
ing  flow  range  and  have  Reynold’s  numbers 
greater  than  unity.  The  low  value  of  D3  = 
1.26  may  result  from  exceeding  the  length 
scales  applicable  to  smaller  aggregates.  This 
lower  value  also  may  be  related  to  the  fact 
that  aggregates  examined  by  Alldredge  and 
Gotschalk  were  examined  in  situ  and  there¬ 
fore  were  not  disrupted  by  handling.  The 
value  of  D3  =  1.26  ±0.06  (±SD)  obtained 
from  the  settling  velocity  data  of  Alldredge 
and  Gotschalk  is  slightly  lower  than  the  val¬ 
ue  of  D3  ■*  1 .39±0. 1 5  (±SD)  calculated  from 
their  gravimetric  analysis.  This  finding  sug¬ 
gests  that  settling  velocity  correlations  may 
underestimate  D3  values. 

Most  of  the  D3  values  calculated  for  ma¬ 
rine  snow  are  smaller  than  the  range  of  val¬ 
ues  observed  for  inorganic  aggregates  formed 
through  diffusion-limited  cluster-cluster  co¬ 
agulation  (DLCA).  There  have  been  pre¬ 
vious  reports  that  biological  molecules  ex¬ 
hibit  different  coagulation  behavior  than 
inorganic  colloids  (Feder  et  al.  1 984).  At  the 
present  time,  we  have  no  specific  explana¬ 
tion  for  these  results.  It  is  generally  accepted 
that  aggregate  characteristics  are  altered  by 
fluid  mixing  intensity,  aggregate  age,  and 
other  factors  not  considered  in  our  analysis 
(Hunt  1986).  Computer  simulation  by 
Meakin  (1988)  demonstrated  that  aggregate 
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restructuring  increases  the  fractal  dimen¬ 
sion.  In  Meakin’s  study,  Z)3  increased  from 
l  .89  to  2.13  when  long  chains  of  an  aggre¬ 
gate  contacted  and  attached  to  other  ten¬ 
uous  structures  within  the  aggregate.  There¬ 
fore,  fluid  mixing  and  turbulence  should 
increase  the  fractal  dimension.  Based  on  our 
current  understanding  of  the  effect  of  fluid 
environment  and  coagulation  kinetics,  and 
Meakin’s  (1988)  results,  we  would  expect 
aggregates  in  the  ocean  to  have  higher  frac¬ 
tal  dimensions  than  those  measured  for 
cluster-cluster,  diffusion-limited  coagula¬ 
tion  and  not  the  lower  values  calculated  in 
this  study. 

We  have  conducted  some  preliminary  ex¬ 
periments  on  the  effect  of  fluid  environment 
on  the  fractal  dimension  of  microbial  ag¬ 
gregates.  When  we  cultured  Z.  ramigera  in 
a  laboratory  reactor  (Omni-culture  reactor, 
Virus  Co.)  aerated  at  1  liter  of  air  min'1  and 
mixed  at  1 60  rpm,  the  fractal  dimension  of 
these  aggregates  was  3.0±0.4  (±SD)  (Wil¬ 
kinson  1989).  This  fractal  dimension  is  es¬ 
sentially  equal  to  the  Euclidean  thiee -di¬ 
mensional  value  for  a  spherical  object.  This 
higher  fractal  dimension  may  be  a  result  of 
more  aggregate  restructuring  in  a  highly 
mixed  vessel  than  occurs  in  less  turbulent 
environments,  such  as  the  ocean.  To  date, 
our  experiments  have  not  had  well-defined 
hydrodynamic  environments  primarily  be¬ 
cause  the  environments  we  are  modeling 
(bioreactors  and  oceans)  are  poorly  char¬ 
acterized  with  respect  to  fluid  mixing  and 
turbulence.  Our  preliminary  findings,  how¬ 
ever,  suggest  that  the  fractal  dimension  of 
aggregates  depends  on  the  fluid  mechanical 
environment. 

We  conclude  that  pure  and  mixed  cul¬ 
tures  of  microorganisms  have  three-dimen¬ 
sional  scaling  properties  typical  of  fractal 
structures.  The  magnitude  of  fractal  dimen¬ 
sions  of  marine  snow  aggregates  is  lower 
than  the  range  expected  for  aggregates 
formed  by  Brownian  motion  of  colloidal  ag¬ 
gregates  and  may  be  related  to  aggregate  re¬ 
structuring  or  aggregate  formation  through 
other  processes,  including  shear  coagulation 
and  differential  sedimentation.  Although 
understanding  of  precisely  how  these  factors 
influence  aggregate  characteristics  is  incom¬ 
plete,  it  is  hoped  that  fractal  geometry  will 


provide  a  basis  for  studying  factors  that  af¬ 
fect  the  structure  of  biological  aggregates  and 
aggregation  kinetics  in  natural  systems. 

Bruce  E.  Logan 
Daniel  B.  Wilkinson 

Environmental  Engineering  Program 
Department  of  Civil  Engineering 
University  of  Arizona 
Tucson  85721 
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